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INTRODUCTION 


ICASE* is operated at the Langley Research Center (LaRC) of NASA by the Universities Space Research 
Association (USRA) under a contract with the Center. USRA is a nonprofit consortium of major U.S. colleges 
and universities. 

The Institute conducts unclassified basic research in applied mathematics, numerical analysis and algo- 
rithm development, computer science, fluid mechanics, and structures and materials in order to extend and 
improve problem-solving capabilities in science and engineering, particularly in the areas of aeronautics and 
space research. 

ICASE has a small permanent staff. Research is conducted primarily by its permanent staff and visiting 
scientists from universities and industry who have resident appointments for limited periods of time as well 
as by visiting and resident consultants. Members of NASA’s research staff may also be residents at ICASE 
for limited periods. 

The major categories of the current ICASE research program are: 

• Applied and numerical mathematics, including multidisciplinary design optimization; 

• Applied computer science: system software, systems engineering, and parallel algorithms; 

• Theoretical, computational, and experimental research in fluid mechanics in selected areas of interest 
to LaRC, such as transition, turbulence, flow control, and acoustics; and 

• Theoretical, computational, and experimental research in structures and material sciences with 
emphasis on smart materials and nanotechnologies. 

ICASE reports are primarily considered to be preprints of manuscripts that have been submitted to 
appropriate research journals or that are to appear in conference proceedings. A list of these reports for the 
period October 1, 2001 through March 31, 2002 is given in the Reports and Abstracts section, which follows 
a brief description of the research in progress. 


* ICASE is operated at NASA Langley Research Center, Hampton, VA , under the National Aeronautics and Space Adminis- 
tration, NASA Contract No. NAS1-97046. Financial support was provided by NASA Contract Nos. NAS1-97046, NAS 1-19480, 
NAS 1-18605, NAS 1-1 8107, NAS1-17070, NAS1-17130, XAS1-15810, NAS1-16394, NAS1-14101, and NAS1-14472. 




RESEARCH IN PROGRESS 


APPLIED AND NUMERICAL MATHEMATICS 


Numerical Simulation of Embedded Streamwise Vortices in a Turbulent Boundary Layer 

Brian Allan 

The use of low-profile vortex generating vanes, which are only a fraction of the boundary layer thickness, 
has been shown to be more efficient than standard-sized vanes, for some passive flow control applications. 
Like the standard size vortex-generating vanes, which are typically the same height of the boundary layer 
thickness, the low-profile vanes still provide effective momentum transfer while having a reduced residual 
drag. These low-profile vanes are currently being investigated for improving engine face flow distortion for 
compact S-shaped inlets. Unlike the standard-size vanes, the low-profile vanes produce streamwise vortices 
which are embedded in the boundary layer. The design of inlets and other aerospace applications, which 
use low-profile vortex generators, need to rely on the ability to model and simulate these devices. This 
investigation will evaluate the ability of a Reynolds averaged Navier-Stokes (RANS) code to simulate an 
embedded stream' wise vortex, generated by a low-profile vane, fixed on a flat plate by comparing numerical 
simulations to experimental data. 

Using a compressible RANS code, developed at NASA, numerical simulations of a single low-profile vortex 
generator vane on a flat plate, have been performed. These simulations were conducted in coordination with 
experiments at NASA Langley’s 20x28-inch Shear Flow Tunnel. Comparisons of the numerical results to 
the experimental data for the single vane, have shown that the turbulence modeling has a significant effect 
on the computed flow field. A comparison between numerical simulations using a Spalart-Allmaras (SA) 
and a Mentor SST model was made to the experimental results. This comparison showed the embedded 
streamwise vortex in the simulations to be more dissipated than was observed in the experiments. However, 
the simulation using the SST model was able to predict the trajectory of the vortex core and the peak 
vorticity, better than the SA model. A comparison of circulation showed that the numerical simulations 
using the SST and SA turbulence models, generated vortices of equal strength. Therefore, the lower peak 
vorticity using the SA model is related to the larger turbulent eddy viscosity predicted by the SA model. 
In general the SST turbulence model was able to predict the trajectory and circulation strength of the 
streamwise vortex. However, both numerical simulations were shown to diffuse the vortex much faster than 
observed in experiments. 

Future work will include detailed experiments using both low-profile and standard-size vanes at different 
device angles on a flat plate. These wind tunnel experiments will be compared to numerical simulations 
providing better insight into the physics and modeling of vortex generating vanes. 

This work was done in collaboration with Pieter Buning, Chung-Sheng Yao, and John Lin (NASA 
Langley). Numerical results were computed using the ICASE PC cluster, Coral. 

Feedback Flow Control Design for Inlet Flows 

Brian Allan 

In order to increase the performance of inlet flows for varying flow conditions, active flow control devices 
for a S-duct inlet design are being investigated at the NASA Langley Research Center, under the UEET 
program. One of the active flow control devices being considered are vortex generating jets. Unlike fixed 
vanes, these devices can be turned on and off when needed. By incorporating the jets with a feedback 
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controller, it is desired that the performance of a S-shaped inlet will be improved over a wider range of 
flow conditions. By using computational fluid dynamics (CFD) the performance of the jets, as well as the 
placement of sensors and actuators, can be explored. CFD will also be used to model, design, and evaluate 
the closed-loop flow control system for the S-duet inlet. 

The inlet flow is simulated using a compressible Reynolds averaged Navier-Stokes code developed at 
NASA. Numerical simulations for the baseline flow of a rectangular duct with an expanding cross section 
have already been performed and compared to experiments conducted at NASA Langley. The numerical 
simulations were able to capture the characteristics of the flow separation pattern but were unable to get 
the exact location of the vortex lift off locations and stagnation points. However, the CFD simulations were 
able to match the surface pressure on the upper and lower surface of the duct for the baseline case. Time 
accurate simulations of vortex generating jets on a periodic ramp have also been performed. These time 
accurate simulations have been used to develop a model of this flow system which was then used in the 
feedback controller design. The objective of this control design was to improve the surface pressure gradient 
on the rainp. Once the feedback controller was designed, it was then simulated by coupling the CFD with 
the control system. The results and experiences from the periodic ramp problem are being used to model and 
design a flow control system for the S-duct inlet which will provide guidance for a flow control experiment 
at NASA Langley. 

Future work will involve baseline simulations of the inlet which wall then be followed by adding flow 
control actuators. Once these simulations are performed the flow system can be modeled and a feedback 
controller designed. These simulations are being performed in conjunction with wind tunnel experiments at 
NASA Langley and will be compared to the wind tunnel results. 

This work was done in collaboration with Pieter Buning and Susan Gorton (NASA Langley). Numerical 
results are computed using the ICASE PC cluster, Coral. 

Real Time Computational Algorithms for Eddy-current-based Damage Detection 

H. T. Banks 

In the held of nondestructive evaluation, new and improved techniques are constantly being sought to 
facilitate the detection of hidden corrosion and flaws in structures such as aeroplanes and pipelines. In 
recent efforts, we explored the feasibility of detecting such damage by application of an eddv-current-based 
technique coupled with reduced order modeling. 

We first, developed a model for a specific eddy current method in which we make some simplifying 
assumptions reducing the three-dimensional problem to a two-dimensional problem (we do this for proof 
of concept). Theoretical results are then obtained that establish the existence and uniqueness of solutions 
as well as continuous dependence of the solutions on the parameters which represent the damage. A least 
squares parameter estimation problem was used in identifying the geometry of the damage. 

To solve the identification problem, an optimization algorithm was employed which requires solving the 
forward problem numerous times. To implement these methods in a practical setting, the forward algorithm 
must, be solved with extremely fast and accurate solution methods. In constructing these computational 
methods, we employed reduced order proper orthogonal decomposition (POD) techniques. This approach 
permits one to create a set of basis elements spanning a data set consisting of either numerical simula- 
tions or experimental data. We developed two different algorithms for forming the POD approximations, a 
POD/Galerkin technique, and a POD/interpolation technique. 

Finally, results of the inverse problem associated with damage detection were obtained using both sim- 
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ulated data with relative noise added as well as experimental data obtained using a giant magnetoresistive 
sensor. The experimental results are based on successfully using experimental data to form the POD basis 
elements (instead of numerical simulations), thus illustrating the effectiveness of this method on a wide range 1 
of applications. In both instances the methods are found to be efficient and robust. Moreover, the methods 
were fast; our findings demonstrate a significant reduction (by factors of order 10 3 ) in computational time. 

This research represents continuing joint effort with Michele L. Joyner (North Carolina State University) 
and Buzz Wincheski and William P. Winfree (NDES, NASA Langley). 

Feedback Control of Nonlinear Systems 

Scott C. Beeler 

We are studying the feedback control of nonlinear systems of ODEs, in particular using the state- 
dependent Riccati equation (SDRE) method. We are especially interested in this technique because it is 
based on optimal control conditions and so its control performance should be very good compared to other 
methods. It is based on the simple form of linear Riccati control solutions while seeking to incorporate 1 as 
much of the nonlinear nature of the problem as possible. This may prove useful in aircraft guidance and 
control problems including control actuator impairment situations. We also are studying the flight dynamics 
of small remote-controlled gliders in an ambient wind field. With unpowered flight the aircraft has a smaller 
velocity, and so the winds have a greater effect on the aircraft's flight properties. This can work in its favor, 
since with efficient, control laws the aircraft can extract energy from thermals or wind gradients to stay 
aloft for long periods of time without needing fuel to run an engine. The research objective is to develop 
new nonlinear control techniques with good performance arid robustness properties, and to study the flight 
properties of small gliders to evaluate their potential for long-term flight. 

We have continued to test our Matlab SDRE feedback control algorithms on example problems with 
state and control nonlinearities, obtaining good results while modifying the several versions of the SDRE 
method we have implemented. An application to a highly nonlinear . inverted pendulum control problem 
is presently being performed, and it has necessitated some alterations in the algorithm to counter the 
more complex nature of this example problem. The glider model has been fine-tuned, and its Matlab and 
Fortran implementations tested. Some very good preliminary results have been achieved for the scenario 
of minimizing altitude loss while traveling from point to point through a wind field. The results achieved 
for both the SDRE and glider codes show the efficiency of the techniques used, and are encouraging as we 
move on to further complex application problems and look to combine these newly developed algorithms 
with established techniques to improve their desirable properties. 

Future efforts include application of the SDRE method to other realistic examples such as aircraft 
guidance problems, possibly a version of the small glider model. It may also be combined with other control 
methods to improve its robustness while retaining its desirable control efficiency. The refined glider model 
will be further tested, using optimal control techniques and later feedback control. Once these small-scale 
control trials have been completed, larger-scale glider surveying problems can be studied, possibly combining 
the small-scale feedback controls with fuzzy logic control for decision-making processes. 

This work is done in collaboration with D. Moerder arid D. Cox (NASA Langley, Guidance and Control 
Branch). 
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Study of the Flux Retrieval Error Behavior with CERES/TRMM Data 

A. Bodas-Salce.do 

The Earth Radiation Budget is a fundamental quantity in climate studies, which on a global scale 
can only be measured with sufficient accuracy from satellites. Radiometers onboard satellites carry out 
directional measurements (radiances), but the physical quantity of interest to obtain the radiation budget 
is the flux density, which is the hemispherical integral of the radiance field. Therefore, the development 
of inversion algorithms (anisotropic models) to retrieve the flux estimate from radiance measurements is 
needed. The purpose of this work is to study the dependency of the error of flux retrieval algorithms with 
respect to the observing geometry. 

In order to carry out this study, following tasks have been developed: 

1. Analysis of 14 days of CERES/TRMM data (CERES stands for Clouds and the Earth’s Radiant 
Energy System, and TRMM for Tropical Rainfall Measuring Mission), by sorting radiance mea- 
surements in angular and cloud cover bins, both in shortwave and longwave bands, and building a 
database at different time scales. 

2. With the previous database, sampling problems that will impact the construction of anisotropic 
models have been investigated. 

3. Shortwave (SW) and longwave (LW) anisotropic models have been built, and a methodology to 
study the performance of the anisotropic models has been developed. 

4. That methodology has been applied to study the flux retrieval error dependency with the observing 
geometry. 

Following this approach, we have demonstrated that a privileged region in the flux retrieval algorithm 
appears in the region of viewing zenith angles around 55 degrees, giving the possibility to obtain top of 
the atmosphere flux estimates with good accuracy and simple anisotropic models by selecting properly the 
observing geometry. 

Although only a few days of CERES/TRMM data have been processed, it seems promising to continue 
this line of investigations by using the whole CERES/TRMM record, in order to minimize sampling problems 
which still may be affecting our results. Furthermore, cloud optical depth should be introduced in the scene 
definition so that a more accurate error analysis may be carried out. 

This research was conducted in collaboration with G.L. Smith (Virginia Polytechnic Institute and State 
University) and E. Lopez-Baeza (Universidad de Valencia, Spain). 

Textbook Multigrid Efficiency for CFD Simulations 

Boris Diskin 

An efficiency goal for the new generation of multigrid flow solvers is to arrive at solutions of the governing 
system of equations in a total computational work that is a small (less than 10) multiple of the operation 
count in one target-grid residual evaluation. Such efficiency is defined as textbook multigrid efficiency (TME). 
TME for the Navier-Stokes simulations can be achieved if different factors contributing to the system could 
be separated and treated optimally, e.g., by multigrid for elliptic factors and by downstream marching for 
hyperbolic factors. An efficient way to separate the factors is the distributed relaxation approach. Design 
of a distributed relaxation scheme can be significantly simplified if the target discretization is factorizable, 
i.e., possessing the following two properties: 

(1) The determinant of the (principally) linearized discrete operator can be represented as a product of 
discrete scalar factors, each of them approximating a corresponding factor of the determinant of the 
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differential Navier-Stokes equations. 

(2) The obtained scalar factor discretizations are stable, easily solvable, and reflect the physical anisotropies. 
Recently, a family of factorizable discretization schemes for the compressible flow equations has been de- 
veloped. The goal of my research performed in the last six months was to analyze these new schemes. In 
particular, I have been focused on discretization of boundary conditions and correct combination of the 
different parts of a multigrid solver to provide a TME solution for the Navier-Stokes equations in different 
flow regimes. 

A detailed discrete analysis has been employed to evaluate several sets of boundary conditions for 
factorizable schemes corresponding to the steady-state compressible Euler equations. Requirements for 
practical boundary conditions have been formulated. Several sets of practical boundary conditions have been 
tested and compared with overspecified boundary conditions. The effects of boundary-condition equations 
on stability and accuracy of the discrete solutions have been analyzed. Explicit correspondence between 
solutions and boundary conditions has been established through a boundary-condition-sensitivity (DCS) 
matrix. Analysis of coefficients of the BCS matrix provides reliable predictions for stability and accuracy of 
the discrete solutions. The following new findings are reported: 

(1) Examples of stable discrete problems contradicting a wide-spread belief that employment of a one- 
order-lower approximation scheme in an O(h )- small region does not affect the overall accuracy order 
of the solution have been found and explained. Such counterexamples can only be constructed for 
systems of differential equations. For scalar equations, the conventional wisdom is correct. 

(2) A negative effect of overspecified (although, exact) boundary conditions on accuracy and stability 
of the solution has been observed and explained. 

(3) Sets of practical boundary conditions for the new factorizable schemes have been formulated. These 
schemes belong to a family of second-order schemes requiring second-order accuracy for some 
numerical-closure boundary conditions. 

The BCS-matrix analysis has also been applied for formulating suitable interface boundary conditions for the 
local relaxation procedure complementing distributed relaxation sweeps in a TME solver for the compressible 
Euler equations. Accurate and stable supplemental boundary conditions for pressure required for solution of 
the incompressible Navier-Stokes equations in the pressure-equation formulation have also been derived by 
means of this analysis. This analysis has been very instrumental for recent TME demonstrations in solving 
collocated-grid discretizations of the Euler and Navier-Stokes equations. 

Currently the efforts are directed to demonstrate TME for general flow computations including general 
boundary conditions and captured shocks. 

This research was conducted in collaboration with J.L. Thomas (NASA Langley). 

Atmosphere Modeling Efforts for 2001 Mars Odyssey Aerobraking 

Alicia Dwyer 

In October 2001, the Odyssey Spacecraft began aerobraking operations at Mars. Aerobraking uses the 
atmosphere of a planet rather than chemical propulsion to reduce the size of the orbit. The success of 
aerobraking depends on the characterization of parameters that affect the safety of the spacecraft. These 
parameters include primarily aerodynamic forces and atmospheric density. Variations in atmospheric den- 
sity constitute the largest orbit-to-orbit uncertainty. Therefore it was desirable to develop models of the 
atmospheric density that accurately predict density at aerobraking altitudes. Predictions of the atmosphere 
influence maneuver decisions and therefore affect the duration of the aerobraking mission and the amount 


5 


of fuel required. 

Throughout the 77 days of aerobraking daily models of the atmosphere were developed. The models were 
based on comparisons between the onboard accelerometer derived density and the Mars Global Reference 
Atmosphere Model (Mars-GRAM) calculated density. Mars-GRAM was also used as the basis for the 
navigation predictions. The density versus time profiles for both the accelerometer and model densities 
were produced for each orbit. Ratios of accelerometer to Mars-GRAM data were obtained for maximum 
and periapsis densities and changes in velocity. Wave models were produced using a least squares fit to a 
selected number of ratios. Because of the large variability in the Mars atmosphere it was not possible to 
rely on a single wave model for extended periods of time (i.e., weeks). For this reason the aerobraking wave 
modeling effort included the analysis of stationary and moving wave models, the effects of using different 
wave numbers, and the optimization of the number of orbits in each fit. The model that produced the best 
fit for a given set of orbits was then used to provide predictions for subsequent orbits. Despite the large 
variability in the Mars atmosphere and the fact that no single model was able to capture the local and daily 
changes for extended periods of time, the modeling efforts proved instrumental in the maneuver decisions 
made throughout aerobraking. 

All the data collected during aerobraking using the accelerometers and the atmospheric modeling efforts 
is being archived and will be used to update current models for future missions to Mars. 

The atmospheric modeling work was done in collaboration with Dr. Robert Tolson and Mr. Gerald 
Keating (George Washington University). 

Development of Aeroassist Capabilities for Planetary Orbit 

Jill Hanna 

Propulsive maneuvers to enter into planetary orbit have been designed since the beginning of space 
exploration. However, to reduce fuel usage and launch mass, more economical approaches for capture into 
planetary orbit have been studied recently. The use of a planetary atmosphere to reduce the orbital energy 
of a spacecraft is the most cost efficient capture strategy to date. One aeroassist strategy, atmospheric 
aerobraking, has been used successfully in three different planetary missions, Magellan at Venus, Mars 
Global Surveyor, and the most recent Mars 2001 Odyssey to change an elliptical orbit to a nearly circular 
orbit. Other methods of aeroassist such as aerocapture are also being studied to reduce the time required 
to modify the spacecraft’s orbital trajectory. 

Mars 2001 Odyssey was the most recent aerobraking success. During the 77-day aerobraking phase 
between October 2001 and January 2002, accelerometer derived density profiles were developed and used 
for atmospheric estimations. The atmosphere was modeled to improve predictive techniques. Three degree 
of freedom (dof) and six dof flight simulations were performed to reproduce and predict the trajectory of 
Odyssey throughout the atmosphere of Mars. Through the six dof work, it was determined that there was 
a significant amount of impingement of the reaction control system thrusters on and around the spacecraft, 
yielding more fuel usage required for the desired effect from the thrusters. The success of Odyssey has had 
a significant impact on future planetary exploration. Odyssey has returned valuable atmospheric data that 
will be useful for future Mars aeroassist opportunities. The science data that is collected from Odyssey will 
aide in choosing landing sites for future missions. The next Mars orbiter is the 2005 Mars Reconnaissance 
Orbiter (MRO). Research is being performed to incorporate the successful Mars Odyssey techniques. The 
research done with MRO will be useful in determining the length of MRO aerobraking arid will improve and 
make more autonomous the aerobraking work supplied by NASA Langley to the Jet Propulsion Laboratory. 
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Other exploration using aeroassist includes a Titan aerocapture mission. The feasibility of a ballute system 
is being studied arid traded with aerocapture using only a heat shield for drag effect . This research will help 
to determine the size, shape, and mass properties of a ballute system. 

Future work includes further improving MRO aerobraking software, changing all software to an updated 
version of programming code, and producing mission runouts and aerobraking statistics. MRO aerobraking 
is scheduled to begin in March 2006 and last for five months. 

Unstructured Nonconforming Multigrid Algorithms for the Solution of Radiation Transport 
Problems 

Kab Seek Kang 

The simulation of radiation transport in the optically thick flux-limited diffusion regime has been iden- 
tified as one of the most time-consuming tasks within large simulation codes. Due to multimaterial complex 
geometry, the radiation transport system must often be solved on unstructured grids. A variety of finite 
volume and conforming finite element schemes have been employed in research at ICASE and at numerous 
national laboratories on this so-called “radiation diffusion” system. Meanwhile, multilevel nonconforming 
finite element or covolume methods have proven flexible and effective on incompressible fluid flow problems. 
The objective of this research is to investigate the behavior and the benefits of unstructured nonconforming 
multigrid algorithms used as linear solvers (inside of Newton approaches), directly as nonlinear solvers, or 
simply as preconditioners in solving steady and unsteady implicit radiation diffusion problems. 

We have completed debugging code 1 of linear versions of a two-dimensional unstructured nonconforming 
multigrid algorithm and have begun to write a report about the behavior and the benefits of the unstructured 
Pi nonconforming finite element method and comparisons with some linearization methods and paralleliza- 
tion. Wo use the P\ nonconforming finite element space on triangular mesh and nested-mesh subdivision to 
automatically generate a. sequence of unstructured meshes on the domain and found the Pi nonconforming 
finite element method of radiation transport problem on inhomogeneous domains resolve very sharp changes 
of energies near the layer where the two different materials meet. 

We plan to measure its scalability on large numbers of processors on the ICASE Beowulf system and 
other machines, and extend into three-dimensional settings. 

Robust Optimization Under Uncertain Design Conditions 

Wu Li 

The ultimate goal is to develop next-generation design tools that are insensitive to variability or un- 
certainty in design specifications. One important example is aerodynamic shape optimization for achieving 
drag reduction under uncertain flight conditions. Our objective is to use both deterministic and probablis- 
tic strategies in optimization process to avoid undesirable decisions made by an optimizer due to a lack of 
information under the uncertain flight conditions. 

A combination of the expected value optimization and the profile optimization is developed to give 
users trade-off options between a greedy reduction of the average drag or a consistent reduction of the 
drag over the range of the uncertainty parameter. This new method is called the constrained multipoint 
approximation method (CMAM) for robust airfoil optimization. CMAM is tested on a lift-constrained drag 
minimization problem for a two-dimensional airfoil in Euler flow, which is formulated with 20 free design 
variables. Each free design variable is the ^-coordinate of a control point in the B-spline representation of 
the airfoil shape. An unstructured grid computational fluid dynamics code, FUN2D, is used to compute 
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the lift/drag coefficients arid their gradients with respect to airfoil shapes and angles-of-attack. Hermite 
spline interpolation formula is used to estimate the expected value of the drag over the Mach range. A 
parallel implementation of CMAM using MPI library reduces the walltime of each optimization iteration to 
one function and one gradient evaluation no matter how many design points are used. Special modifications 
of the standard trust region optimization strategy are used to reduce potential random shape distortions 
during the optimization process. Also, due to difficulty or impossibility of obtaining accurate derivatives 
for aerodynamic quantities, we are also working on a derivative- free optimization (DFO) method for robust 
airfoil shape optimization. 

In the last two years, researchers at ICASE have developed the expected value optimization method and 
the profile optimization method for robust airfoil optimization under uncertain flight conditions. The ex- 
pected value optimization method uses a dynamic numerical approximation of the expected value of the drag 
to avoid off-design degradation problems encountered by the traditional multipoint optimization method, 
while the profile optimization method uses a smart descent direction to achieve a consistent drag reduction 
over the entire Mach range and search for a robust optimal solution. By selecting various ratios of two input 
control parameters, CMAM can generate optimal solutions corresponding to the multipoint optimization 
method, the expected value optimization method, the profile optimization method, and a combination of all 
these methods, respectively. This shows how these different optimization methods are related to each other. 
The numerical simulation results indicate that CMAM is flexible for testing and comparing various robust 
optimization strategies. By using a special trust region strategy, it generates intermediate iterates without 
any undesirable noises in airfoil shapes, no matter what optimization strategy is used. With only four design 
points, CMAM generates smooth optimal airfoils with smooth drag profiles over the given Mach range. Op- 
timization techniques can automate many design tasks, providing better designs in a shorter time. However, 
traditional optimization methods require precise specification of the objectives and constraints. By contrast, 
robust optimization methods allow uncertainty in design specifications. However, when there is a nonlinear 
uncertainty parameter in the optimization model, one could only solve an approximation model of the original 
optimization problem. Each approximation method has its own strength and weakness. Some methods allow 
a designer to explore the design space, while others generate improved designs using the smallest possible 
change in the baseline configuration. It is critical to be able to analyze and compare various techniques for 
robust optimization involving highly nonlinear objective/constraint functions with respect to the uncertainty 
parameter (such as the drag with respect to the Mach number). CMAM uses a general optimization model 
to encompass all existing robust optimization methods for airfoil shape optimization under uncertain flight 
conditions and allows a systematic study of various robust airfoil optimization methodologies. 

CMAM will be extended to study airfoil shape optimization with multiple uncertainty parameters. In 
collaboration with ICASE, the mathematical properties of the method will be analyzed in order to generalize 
the method for use in future MDO applications. In collaboration with Michael Wagner, we will also develop 
DFO methods to deal with more realistic and complex aerodynamic shape optimization problems. 

This research was conducted in collaboration with Sharon Padula (NASA Langley) and Michael Wagner 
(Old Dominion University and ICASE). 

The Coral Project 

Josip Loncaric and Manuel D. Salas 

The cost of developing complex computer components such as CPUs has become so high that scientific 
applications alone cannot carry the full burden. Scientific computing needs to use mass market leverage 
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to overcome the cost barrier. A cost-effective alternative to high-end supercornputing was pioneered by 
Beowulf, a duster of commodity PCs. By now, very high performance Beowulf clusters can be built using 
fast commodity PCs and switched Fast Ethernet. We are exploring the benefits and the limitations of this 
approach, based on applications of interest to ICASE. 

The initial phase of the Coral project, consisting of 32 Pentium II 400 MHz nodes and a dual-CPU 
server, demonstrated aggregate peak performance in excess of 10 Gflop/s, with sustained performance on 
CFD applications of about 1.5 Gflop/s. In order to provide a richer environment for further experimentation, 
a dual-CPU configuration was chosen for the second phase of the Coral project. We have added 1G dual 
Pentium III 500 MHz machines and two dual-CPU file servers. The third stage of this project added 16 dual 
Pentium III 800 MHz machines and a 32-node low latency 1.25 Gbps Giganet cLAN network fabric. This 
year the project is in its fourth stage, and 24 of the original 400 MHz machines were replaced with much faster 
1.7 GHz machines. The resulting system contains 96 compute CPUs and 7 server CPUs with an aggregate 
of nearly 55 GB of RAM and 2 TB of raw disk space. Coral has an excellent price/performance ratio, almost 
an order of magnitude better than an equivalent proprietary supercomputer design. This conclusion is based 
on our experience with a variety of applications, ranging from coarse-grained domain decomposition codes 
to communication- intensive parallel Tenderers. Coral has computed solutions to many large scale problems, 
and users needed a filesystem large enough to store their output. Therefore, we have converted one of the 
original 400 MHz machines into a file server with a RAID-5 array consisting of four 160 GB ATA/133 disks. 
The formatted capacity of this new filesystem is 481 GB (where GB=10 9 bytes). This RAID-5 array can 
read data at 66 MB/s. The file server has been observed to deliver up to 35 MB/s when serving multiple 
nodes via the Gigabit Ethernet interface using the NFS protocol. The incremental cost of this file server was 
only S2.24/GB of filesystem space. 

We will continue to use this cluster to develop and run research codes of interest to ICASE and NASA 
Langley, and to evaluate price/performance tradeoffs among various hardware, software, and networking 
configurations. 

Active Shielding and Optimal Control of Environmental Noise 

Jo sip Loncaric and Semyon V. Tsynkov 

Rejection of exterior noise caused by periodic sources such as propellers or turbines would significantly 
enhance passenger comfort and reduce noise fatigue on long flights. Passive sound absorbing materials help at 
high frequencies, but, to be effective below about 1 kHz their weight penalty becomes significant. Active noise 
control can reduce low frequency noise with less weight penalty. Based on the mathematical foundations of a 
new active technique for control of the time-harmonic acoustic disturbances, we have developed a numerical 
technique which can suggest good locations for sensors and actuators. 

Unlike many existing methodologies, the new approach provides for the exact volumetric cancellation 
of the unwanted noise in a given predetermined region of space while leaving those components of the total 
sound field deemed as friendly unaltered in the same region. Besides, the analysis allows us to conclude that 
to eliminate the unwanted component of the acoustic field in a given area, one needs to know relatively little: 
only the perimeter data (the total acoustic field and its normal derivative) are required. The mathematical 
apparatus used for deriving the general solution is closely connected to the concepts of generalized potentials 
and boundary projections of Calderon's type. This exact general solution can be computed at polynomial 
cost, and good actuator locations determined via a procedure which progressively restricts the locations where 
actuators may be placed. The answer depends on the chosen optimality criterion. While L 2 optimality of 
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control inputs is typically sought, this criterion is not physically meaningful. Minimizing the L\ norm of the 
control effort has a clear physical interpretation, but it is very difficult to solve numerically. Fortunately, our 
numerical results for the two-dimensional case and our analytic proof for the 1-D case show that perimeter 
controls are L\ optimal. This is a strong basis for our conjecture that perimeter controls are L\ optimal 
in general. Finally, we've shown that noise control can absorb power. The power optimal controls load the 
noise sources to increase* their power output, then absorb half of the power increase as the other half radiates 
to infinity. 

These results on L 2 optimality, L\ optimality and power optimality provide useful guidance for the exact 
volumetric cancellation of noise. We plan to seek funding for further work on approximate noise cancellation, 
as well as other possible applications which may include different physics, such as electrodynamics, and 
different formulations of the boundary-value problems, such as scattering. 

It is our pleasure to acknowledge most useful discussions with Jan Hesthaven, Wu Li, Michael Overtoil 
and Victor Rvanben'kii. 

Human/Robotic Exploration of the Solar System 

Josip Loncartc and Manuel Salas 

Exploration of the solar system will be most effective if both humans and robots are synergist ically 
combined. Done correctly, this approach can reduce risks, improve efficiency, and accomplish goals faster. 
The challenge is to understand the ways in which this could be accomplished and how this mix might evolve 
over the next 10 40 years with the incorporation of revolutionary aerospace systems concepts. We have 
organized a workshop which brought together experts from academia, government, and industry to address 
these challenges. 

This workshop, held on November 6 7, 2001, was aimed to support the development of a preliminary 
plan which would maximize the scientific return and enable the human exploration and development of space. 
The scope of this effort includes both planetary science and “in space 1 ' platform science applications beyond 
low Earth orbit. Specific objectives are: (1) identification of advanced revolutionary systems concepts, 
(2) identification of required technologies to enable these capabilities, (3) an evaluation of the evolution 
of the relative roles of humans and machines to implement these concepts, and (4) an identification of 
the science that would be enabled by these capabilities. Support for this project was provided by the 
Revolutionary Aerospace Systems Concepts (RASC) activity at NASA Langley. The workshop attracted 
93 participants from government, academia and industry. After the workshop, session leaders gathered at 
ICASE on November 8th, 2001 to draft the final report. Moreover, the presentations and the NIAC RFI 
submissions were converted to PDF format and published on the workshop web site. 

Additional visual illustrations of the exploration concepts presented at the workshop will be prepared 
by JF&A. The final NASA Conference Proceedings should be published in 2002. 

This work is being conducted in collaboration with Robert Cassanova (NIAC) and Lewis Peach (USRA). 

Aerodynamic Transonic Drag Prediction Using an Unstructured Multigrid Navier-Stokes Solver 

Dimitri J. Mavriplis 

The objective of this work is to assess the capability of an unstructured grid Navier-Stokes solver at 
predicting drag in the transonic cruise regime of a transport aircraft. 

The flow over a wing-body transport aircraft configuration was computed on a grid of 1.6 million points 
and a grid of 13 million points at various Mach numbers and lift coefficient values. The results were compiled 
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as a series of drag polars and drag rise (drag versus Mach number at fixed lift coefficient.) curves. A total of 
72 individual cases were computed on the 1.6 million point grid, which was performed in about one week on 
the ICASE PC cluster using 16 to 32 processors. A total of six cases were computed on the 13 million point 
grid which was performed on the SGI NAS Origin 2000 machine using 64 to 128 processors. The drag results 
were found to be in relatively close agreement with results from other validated structured, unstructured, 
and overset grid method codes in the workshop, although a consistent bias in the lift at a given angle of 
attack existed between these numerical results and the experimental wind tunnel values. The ability to 
compute large numbers of cases with good accuracy on unstructured grids using inexpensive PC clusters was 
demonstrated through this exercise. Results from this study were presented at the 40th AIAA Aerospace 
Sciences Meeting in Reno, NV, January 2002. 

More computations on additional grids are to be performed and the ability of adaptively generated 
meshes to reduce solution time and increase accuracy are also to be studied. 

Investigation of Multigrid-implicit Higher-order Accurate Time-stepping Procedures for 
Unstructured Mesh Navier-Stokes Simulations 

Dimitri J. Mavriplis and Giridhar Jothiprasad 

The accurate solution of transient fluid flow problems as required in large eddy simulations, can become 
very cpu-time intensive due to the long simulation times required and large number of time steps employed. 
This work explores the feasibility of reducing the overall simulation time for a given level of accuracy by 
resorting to a smaller number of higher-order accurate time steps and using more efficient implicit solution 
techniques at each time step. 

Based on the work of H. Bijl, M. Carpenter, and V. Vatsa, an implicit Runge-Kutta time stepping 
scheme has been implemented in an existing two-dimensional unstructured Navier-Stokes solver. A four- 
stage and a six-stage Runge-Kutta scheme are employed which are third-order and fourth-order accurate in 
time, respectively. Each stage of the Runge-Kutta schemes represents an implicit nonlinear problem which 
is solved using the agglomeration multigrid scheme previously developed for steady-state problems. This 
multigrid solver has also been used as a solver for a second-order backwards difference time stepping scheme. 
Design accuracy for both Runge-Kutta schemes has been demonstrated for laminar flow over a circular 
cylinder. A comparison of the time required to solve the nonlinear problem at each time step using a linear 
multigrid algorithm and a nonlinear multigrid algorithm has also been investigated. The linear multigrid 
algorithm has been found to be substantially faster than the nonlinear multigrid algorithm for this task. 

The relative efficiencies for a given accuracy level of the various Runge-Kutta schemes and backwards 
differencing schemes solved with linear and nonlinear multigrid methods are to be compared. The effective- 
ness of using the multigrid solvers as preconditioners to a Newton-Krylov method is also under investigation. 
These results are to be presented at the AIAA Fluid Dynamics conference in June 2002. 

Parameter- varying Gain-scheduled Fault Tolerance Control of an Aircraft 

Jong-Yeob Shin 

Recent aircraft flight control design research has concentrated on the development of control systems 
which stabilize an aircraft robustly and allow to recover aircraft performance under adverse conditions 
(actuator/sensor/engine failures and bad weather). One promising approach to designing robust Fault 
Tolerance Controller (FTC) for an aircraft under adverse conditions is use of linear parameter varying 
(LPV) control synthesis, since aircraft dynamics vary significantly according to angle of attack, velocity, 
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dynamic pressure, and so on. Aircraft dynamics can be described by a quasi-LPV model with scheduling 
parameters which are also states of the dynamics model. 

To apply an LPV control synthesis to design a FTC, control component failure cases should be modeled 
as an LPV system. Usually, control surface failure cases are modeled as linear functions of fault parameters. 
The fault parameters should be measurable on-line to implement an LPV controller. We have currently used 
the optimal two-stage Kalman Filter to estimate fault parameters. In this case, the scheduling parameter is 
estimated with bounded uncertainty. We have formulated an LPV synthesis methodology with scheduling 
parameter uncertainties. The scheduling parameter uncertainties are represented as an uncertainty block in 
the LPV synthesis methodology. The designed LPV controller may lead a conservative result since the LPV 
synthesis methodology can not count on the uncertainty structure. In order to reduce conservatism of the 
controller, scaling factor 5 is introduced. The robust LPV controller problem is to find S and a controller 
K\ which is not a linear matrix inequality problem. We have completed coding the robust LPV controller 
problem in MATLAB using an iteration method. 

In future work, the robust LPV controller will be integrated with the fault parameter estimator. Also, 
a model of actuator failures will be considered as nonlinear functions of the fault parameter. Developing 
as reliable quasi-LPV model of an aircraft dynamics with fault parameters is still an open problem. The 
problem will be investigated based on research of a function substitution method. 

This research is conducted in collaboration with Dr. Christine Belcastro (NASA Langley) and Prof. X. 
Eva Wu (Binghamton State University of New York). 

Knowledge Discovery in Large Scale Distributed Data Sets with High- and Infinite-dimensional 
Data 

Eduardo Socolovsky 

To produce new approaches, algorithms, and tools to help extend the limits of computational feasibility 
and reduce the cost of performing clustering and dimension reduction on high-dimensional large-scale dis- 
tributed data sets. Most presently used data analysis techniques become unfeasible or costly for very modest 
size data sets. The key factors to overcome this problem, are the computational efficiency of the algorithms, 
as well as the scalability of the algorithm and its implementation. Surveys indicate that between 70% and 
80% of all data mining applications involve some form of clustering, which is a critical step in, e.g., statistical 
data analysis, pattern recognition, KDD, image processing, and data layout. Clustering high-dimensional 
data sets is a contemporary challenge. Dimension reduction is often crucial for high-dimensional data since, 
generally, the computational complexity grows exponentially with the dimension, and is also important as 
an analysis technique, such as Principal Component Analysis (PC A), also known as empirical orthogonal 
functions (EOF), Karhunen-Loeve transform, or normal modes. 

In clustering a measure of similarity and/or a measure of density are used to detect groups or patterns in 
data. The most common is the Euclidean distance between data points, and the standard density is “number 
of points per unit volume/ 1 However, for high-dimensional data, the Euclidean distance does not work well, 
and the above density becomes useless since the required number of points grows exponentially with the 
dimension. For higher dimensional data and pattern recognition, clustering methods adopt more adequate 
measures of similarity, like Pearson’s product-moment correlation or the cosine measure. An important 
problem with the above and other similarity measures in high dimensions is that, the triangle inequality 
doesn’t hold. Among the results we have obtained, is a measure of dissimilarity strongly connected to the 
cosine measure that satisfies the triangle inequality. With the correlation similarity measure the bulk of the 
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computation and error propagation is in the dot-products, in particular the O(N) inner products required for 
the similarity measure matrix. We outlined new mixed agglomerative-partitioning algorithms that (at least 
in the average case) compute only O(N) inner products using our results on the two measures to generate 
approximate similarity and dissimilarity matrices (or equivalent hounds), that can be obtained from already 
computed inner products plus a few more floating point operations. 

In PC A the set of variables or features are transformed into a new mutually uncorrelated (orthogonal) set 
of variables, and a smaller representative subset of variables is selected with a variance maximization criterion. 
Computationally, the transformation is obtained from the Singular Value Decomposition (SVD) of the matrix 
formed by the data points, and the new variables are given by the eigen- vectors corresponding to the largest 
singular values of the matrix. The original data is then projected onto the space spanned by the selected 
eigen-vectors. However, finding the SVD is computationally intensive and is a contemporary challenge to 
adapt for distributed data, which is resolved by adopting a PC A approximation. Furthermore, for some data 
sets the new variables are difficult to interpret and PCA does not provide a satisfactory representation, as 
frequently occurs with climate data sets. Often a representation in term of the original variables is preferred, 
hence as is done in Independent Component Analysis, we intend to investigate the use of other directions. 
The new directions are obtained or suggested by the results of clustering the original variables, thus obtaining 
a different matrix factorization. For scalability, current data mining system architectures favor a distributed 
agent approach to generate local components that are later merged into a global result, often based on a 
peer-to-peer agent collaboration and negotiation which can result in high communications demands. We will 
also explore Data Swarming, a non-traditional biologically inspired approach to software design and scientific 
applications in which the results emerge from the work of a collection of agents that perform simple tasks 
and communicate through a global medium. 

Review of Fiision Propulsion for Spacecraft 

Gerald Walberg 

This work was carried out under the Revolutionary Aerospace Systems Concepts (RASC) program. 
Direction was received from NASA Headquarters to include fusion propulsion in the FY 2002 R ASC studies. 
In order to address fusion propulsion in an informed manner, it was decided to carry out a literature survey 
of fusion research and of spacecraft designs using fusion propulsion in order to assess the current state of the 
art. 

The literature surveyed consisted of recent (1980 to present) survey articles, research results, proposals 
for proof-of-concept investigations and published spacecraft designs. Both magnetic-confinement fusion and 
inertial-confinement fusion were addressed. It was found that a vast body of experimental and analytical 
research results exists and is sufficient to allow credible conceptual designs of fusion propulsion systems. 
When spacecraft are designed based on extrapolations of this currently existing database, they are very large 
(typical overall dimensions of 100 200m) and massive (6000-7000 metric tons), and require extensive in-space 
assembly. In order to achieve acceptably small and compact fusion propulsion systems, new approaches to 
fusion that have not yet been experimentally validated must be employed. 

A RASC White Paper has been prepared and will be used to set the acceptable parameter space for the 
2002 studies. 
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Methods Of Separation Assurance For Instrument Approach Procedures 

Maria Consiglio 

A method to provide automated air traffic separation assurance services during approach to or departure 1 
from a non-radar, non-towered airport environment is investigated. The method is constrained by the current 
NAS infrastructure without radical changes or ambitious technology investments. 

The proposed procedures are designed to grant access to a large number of airfields that currently 
have no or very limited access under Instrument Flight Rules (IFR) thus increasing mobility with minimal 
infrastructure investment. While mobility, defined as point to point transportation of people and goods may 
be increased total distance traveled and demand on busy airports may be reduced. 

Following the definition of the procedures a simulation study is under way that models the ground-based 
airport automation to estimate the system behavior and response under given conditions. We completed 
coding the first release of the simulation driver and sequencing algorithm. We are currently doing testing 
and validation and verification of simulation results. 

This work is done in collaboration with Sheila Conway (NASA Langley). 

Analysis of Causal Event Relationships in a Distributed Air Traffic Management Simulation 

Maria Consiglio 

Understanding the behavior of a parallel, distributed system can, among other things, lead to more 
efficient implementations. Distributed systems implementations that assume a global clock are often simpler 
to implement and easier to understand. Potential inconsistencies and deadlock conditions due to the inherent 
non-determinism of a distributed computation are prevented as the computation proceeds time stepping, 
synchronizing at every tick of the global clock. An alternative approach allows distributed components of 
the simulation to use their local clocks while ensuring event ordering as imposed by the causality constraints. 
Conservative and optimistic (time-warp) implementation methods can be used once the system behavior has 
been fully characterized. 

Identifying the causal and temporal relationships can eliminate the need of an artificially imposed global 
time which can introduce unnecessary communication overhead, reduce the degree of parallelism, limit the 
scalability of the system, and compromise the validity of the computation. 

In this work, the causal and temporal event dependencies in a distributed air traffic simulation (FFSim) 
are analyzed to verify the viability of a distributed time implementation that ensures correct results and sat- 
isfies repeatability requirements. FFSim (Free Flight Simulation) is a workstation- based, medium fidelity air 
traffic simulation designed to test DAG-TM (Distributed Air /Ground Traffic Management) operational con- 
cepts that is being developed in the Air Traffic Operations Laboratory (ATOL) at NASA Langley Research 
Center. 

This work is done in collaboration with Sheila Conway (NASA Langley). 

Tidewater Research Grid Project (TRGP) 

Thomas M. Eidson and Josip Loncaric 

A recent trend in scientific computing has been the increased use of Grid computing. Grid computing 
is defined as the development and execution of distributed applications across wide-area networks where ad- 
ministrative and security issues are non- trivial. Very high-performance computers have become increasingly 
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expensive and researchers are looking to the use of large numbers of lower-performance machines to meet 
their computational needs. Unfortunately, these machines are located at various locations where the Inter- 
net (or an internet) is involved, thus the need for Grid computing. A number of projects around the world 
are focusing on the technology and security of remote job execution, few are looking at the security issues 
involved with managing a distributed user base located at multiple sites. TRGP was formed to assist ICASE 
and other related organizations in learning Grid software basic technology and to provide a base platform 
for the development of Grid programming environments. Additionally, the user- management security/trust 
issues are being studied. 

The current focus of the project has had two objectives. The first was to use currently available software 
to implement a Grid infrastructure. The Globus Toolkit w r as used to provide the distributed computing 
functionality. This includes remote process management, file transfer capability, message passing and other 
signaling, arid security. The PGP system w r as used to add secure person-to-person communication since 
the participants in TRGP will be located at separated sites. These packages have been installed and a set 
of Perl scripts were written to provide a simple arid consistent interface for TRGP users. This has been 
all documented via Web pages. The goal being that both system administrators and users wall have single 
dear source to learn about TRGP and Grid computing. The second objective was to recruit organizations 
to part of TRGP. A workshop was held to present the above TRGP system to potential organizations. 
Representatives from Old Dominion University, Jefferson Labs, and W illiam and Mary were invited. The 
response w r as favorable as all have expressed interested in joining TRGP. 

In the next phase of this project we have three objectives. The first will be to give installation support 
to any organization that wishes to join TRGP. While the goal of the TRGP documentation and software is 
to be self-sufficient, it will be important to get feedback from early users to improve the documentation and 
software. The second objective will be to work with users to develop Grid applications. The large number 
of parallel applications designed for homogeneous parallel architectures do not necessarily port to Grid 
architectures. Applications will need to have a heterogeneous nature to map to the heterogeneous collection 
of machines and networks. One w r ay this can be done is to try ncw : algorithms. However, a more fruitful 
approach may be to develop more complex, integrated applications. This can be done b\ coupling data 
analysis to an existing application so that application can be interactively monitored and possibly steered. 
Another approach is to develop applications that solve more complex problems and involve multiple types 
of physics. Multi-dimensional optimization is one example of this type of application. The final approach 
in this phase will be to study the Globus Toolkit, and relating programming techniques. The objective will 
be to become familiar with any new r heterogeneous, distributed programming techniques to provide the best 
assistance to TRGP users. 

15th International Conference on Theorem Proving in Higher Order Logics (TPHOLs 2002) 

Cesar Munoz and Hanne Gottliebsen 

Theorem proving technology is currently being used in industry and government in places such as Intel, 
Rockwell-Collins, Honeywell Aerospace, NASA Langley Research Center, NASA Ames Research Center, 
and the Naval Research Laboratory. Several universities in the United States and around the world also 
have an active formal methods research program to include theorem proving and applications. Industrial 
usage and increase in system complexities are placing new r demands on theorem proving, such as improved 
automation, domain specific theories, and interoperability of multiple verification tools. For technology 
transfer to succeed, theorem proving researchers need to understand these new demands and create new 
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research thrusts to address the issues. 

TPHOLs is the premier conference on the theory, development, and application of higher order logic 
and theorem proving. The conference brings together international experts from government, industry, and 
academia for four days of research presentations, work-in-progress presentations, discussions, and demon- 
strations. In 2001, ICASE, NASA Langley, and Concordia University put together a bid to host the 15th 
edition of this prestigious conference on August 20- 23, 2002 in Hampton. VA. The bid was chosen in a voting 
process by the majority of the TPHOLs community. Since then, the TPHOLs 2002 organizing committee 
has customized and improved the START Conference Manager Web Server to support: electronic paper 
submissions and registrations; and web-based reviewing of submitted works and voting for the selection of 
the TPHOLs 2003 site. 

The 15th edition of TPHOLs will be co-located with the First Workshop on Formalizing Continuous 
Mathematics. More than 60 researchers in the area of theorem proving and higher-order logics are expected 
to attend both meetings. The proceedings of the conference will be published in the Springer- Verlag's Lecture 1 
Notes in Computer Science series. 

The conference and workshop are being organized with the collaboration of Victor Carreno (NASA 
Langley), Sofiene Tahar (Concordia University), and Andrew Adams (University of Reading). 

Fundamental Algorithms for the En-route Phase of DAG-TM 

Cesar Munoz , Alfons Geser , and Hanne Gottliebsen 

Free-Flight is a new concept for air traffic management. Under free-flight rules, air traffic separation is 
distributed among the actors of the aerospace system. DAG-TM is the NASA implementation of the Free- 
Flight concept. It consists of a set of elements , also called Concept Elements , each one of which addresses 
a particular phase of a flight. Concept element 5 corresponds to the en-route phase. The main problem in 
this phase is the avoidance of air traffic conflicts. The objective of this work is to develop a set of formally 
verified algorithms for the en-route phase of DAG-TM. Formal verification will substantially increase safety 
of the increasingly autonomous and complex new air traffic management systems. 

In the past, we have successfully applied formal methods techniques to design a three-dimensional 
algorithm for tactical conflict resolution (KB3D) and strategic conflict resolution and recovery (RR3D). In 
this work, we have extracted and modified the conflict detection algorithm in KB3D towards an algorithm. 
CD3D, that is not only correct but also complete. We have so developed a suite of fundamental algorithms 
for the en-route phase of free flight that includes: 

• CD3D: A correct and complete three-dimensional conflict detection algorithm. 

• RR3D: A strategic algorithm for conflict resolution and recovery that is based on geometry. 

• KB3D: A correct , improved version of the original conflict resolution algorithm. 

The RR3D algorithm has significantly improved and matured. Its current version and a mathematical 
argument for its correctness are documented in an ICASE Technical Report. We established two theorems 
that we believe represent a breakthrough in our attempt to prove correctness of RR3D. All three algorithms 
have been analyzed using rigorous mathematical techniques. Furthermore, CD3D and KB3D have been 
mechanically verified in the PVS system. Prototype implementation of the algorithms are available through 
the Internet and simulations on high-density free-flight scenarios are being performed. 

Formal methods are a promising technology for safety of avionics systems. We will continue the de- 
velopment of fundamental algorithms for other phases of DAG-TM using formal methods techniques. In 
particular, we are beginning to study the merging and self-separation problems of Concept Element 11 
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(Terminal Area). 

This work was done; in collaboration with Gillos Dowek (INRIA) and Florent Kirchncr (ENAC). 

A Formal Verification of the SPIDER Diagnosis Protocol 

Alfon . s Genet 

SPIDER (Scalable Processor-Independent Design for Electromagnetic Resilience) is a case study for the 
Federal Aviation Agency that shall demonstrate the feasibility of Formal Methods for the certification of 
safety-critical hardware according to document DO-254. It shall provide feedback on problem areas and 
shall serve as material for DO-254 training. 

SPIDER is composed of a number of Processing Elements (PEs) that are connected via the Reliable 
Optical Bus (ROBUS). The ROBUS consists of two columns of nodes, the BIUs and the RMUs, which 
are connected to each other but not connected among their own kind. The nodes are fault-containment 
units, i.e., all nodes fail independently. Nodes may exhibit various kinds of faulty behavior, including two- 
faced, malicious behavior. To counteract the danger that emergence of faults confuses the system, each 
node maintains and collects evidence about the health status of other nodes through direct and indirect 
observations. The purpose of the Diagnosis Protocol is to distribute each node’s knowledge about the health 
status reliably, and to form convictions, i.e., correct claims of a node being faulty that all good nodes share. 
Diagnosis works only if enough hardware is good (Maximum Fault Assumption, MFA) and the health status 
information is intact (Correct Active Sources). Having finished the correctness proof in PVS at the start 
of the period, our primary goal was to simplify and to document the proof. Opportunities to migrate to a 
dynamic readmission of transiently failed nodes were discussed. During this discussion a previously neglected 
important point came up: can we prove that diagnosis strictly increases the reliability of the ROBUS? We 
found that for this purpose the MFA has to be weakened such that only undeclared nodes are counted. 
This dynamic MFA does indeed profit from convictions as can easily be shown. The Diagnosis Protocol, the 
Interactive Consistency Protocol, and the correctness proofs of the two needed a few essential changes to 
accomodate the dynamic MFA. We have implemented these design changes, including the correctness proofs, 
and are now about to write an ICASE Technical Report. 

In the near future we plan to start designing a Readmission Protocol that allows formerly faulty but 
recovering nodes to be reconsidered as good nodes. 

This work was done in collaboration with Paul S. Miner (ATM, NASA Langley). 


17 


FLUID MECHANICS 


Seedless Laser Velocimetry with Laser-induced Thermal Acoustics 

Roger Hart 

Noil-intrusive measurement of flow velocity by various laser-based methods provides information of great 
significance to the experimental fluid dynamics community. However, the only methods to have achieved 
widespread use to date are laser doppler velocimetry and particle interval velocimetry, both of which require 
the flow to be seeded with small (I micron) particles to serve as light scatterers. Seeding is not feasible 
in some wind tunnels due to concerns over removal of spent seed, dogging of flow straightening screens, or 
abrasion of finely polished surfaces. Additionally, there are regions in air flows of interest, such as vortex 
cores over delta wings or in recirculation regions behind rearward facing steps or leading-edge slats, where 
useful seed concentrations may be difficult or impossible to achieve. There is, thus, considerable interest in 
a laser velocimetry technique that can provide the dependability, ease of use, and quality of data of existing 
methods without seeding. Laser-induced thermal acoustics (LITA) is a relatively new technique which has 
great promise as a robust, reliable seedless laser velocimetry method. Our objective is to further develop 
the theoretical and technical basis for LITA velocimetry in the laboratory, and to demonstrate the utility of 
LITA velocimetry through a variety of Teal-world' wind tunnel measurements at NASA Langley. 

Considerable theoretical and laboratory evaluation of LITA velocimetry has already been accomplished, 
although work continues in support of the improvement of our designs for fieldable measurement systems. 
Our overall design philosophy springs from the realization that a system must be relatively easy to install and 
use and must reliably produce dependable data; a system that does not meet these requirements will be of no 
interest to our user community. To that end our designs stress modularity, easy transport and installation, 
and the use of compact, turn-key laser systems. Of great significance as regards stability and accuracy is 
our use of a novel grating-based quadrature demodulation technique (invented by us). The current reporting 
period has been devoted to the completion of a two velocity component LITA velocimeter incorporating 
lessons learned from our earlier one-component prototype. The new instrument is fully remotely operable 
from a wind tunnel control room, and should be significantly more stable than the earlier version. 

We are currently scheduled to conduct a demonstration of the new velocimeter at BART beginning 
in early May 2002 in which we will make velocity measurements of the flow above a high-lift wing. Our 
data will be compared with existing PIV data sets for that configuration. In February 2003 we will begin 
our first demonstration of supersonic LITA velocimetry at PCT, investigating the flow around a diamond- 
shaped model at free stream velocities ranging from transonic to Mach 3.5. The possible application of 
LITA velocimetry to the high-pressure exhaust flow from a (leftover) NASP model at 15” Mach 6 is under 
discussion. 

This work is pursued in collaboration with R.J. Balia and G.C. Herring (AMDB, NASA Langley). 

Time Accuracy of Large Eddy Simulation for Unsteady Flows 

Guowei He 

The objective of this research is to develop subgrid scale (SGS) modeling for Large Eddy Simulation 
(LES) of unsteady flows. The previous research has focused on the ability of LES on single-time flow proper- 
ties such as energy spectrum. Attention has never been paid on the two-time or frequency energy spectrum. 
In fact, the frequency spectrum has wide application to unsteady flows. For example, the frequency spectrum 


18 


I 



of velocity fluctuations determines the frequency distribution of far-field acoustic energy in aeroacoustics: 
the frequency content of turbulent energy characterizes the turbulent disturbance environment for receptiv- 
ity; the noise reduction, synthetic jets, and coherent structure in flow control are associated with frequency 
energy spectrum. 

We have shown that the eddy viscosity SGS models over-predict the two-time energy spectrum. The 
over-prediction could In; very significant. Further, the necessary conditions for SGS models to correctly 
predict time correlation are derived. These conditions are used to investigate SGS models and suggest 
that, the random forcing model and Lagrange dynamical model are recommendable for the two-time energy 
spectrum. They also indicate that the coefficients in the eddy viscosity SGS models should be history- 
dependent. We have developed a new two- parameter model and sub-Gaussian model for the two time energy 
spectrum. The former is constructed by the “sweeping” and “oscillation” parameters and the latter by the 
mean velocity of sub-Gaussianality. 

Our future work will develop the Langevin-type and Lagrange dynamical SGS models for time accuracy 
of LES in unsteady flows. 

This work is performed in collaboration with R. Rubinstein (NASA Langley). 

Uncertainty in Computational Fluid Dynamics 

GuoweA He 

The objective of this research is to develop a statistical approach to uncertainty in Computational Fluid 
Dynamics (CFD). We focus on the errors of CFD induced by either turbulence modeling or uncertainty 
in physical models. The motivation of this research is that NASA needs credible CFD simulations and 
aerodynamical optimization needs the quantification of uncertainty in CFD as input parameters. 

We are developing the mapping closure approximation to quantify the uncertainty propagation. The 
mapping relates the random inputs and outputs in deterministic fashion. Actually, the truncated mapping 
has been used to calculate the first and second order moments in sensitive analysis. The demonstration cases 
are the passive scalars in shear flow and reactive scalars in random Gaussian flows. The probability density 
function, conditional dissipation rate, and conditional Laplacian are calculated and are in good agreement 
with direct numerical simulation. It is shown that the approach is able to describe the propagation of 
uncertainty in physical models. Two ICASE reports on this topic have been finished. 

Our future work will apply the mapping closure approximation approach to the quasi-one-dimensional 

Euler CFD. 

Multiple-relaxation-time Lattice Boltzmann Models in Three Dimensions 

Li-Shi Luo 

The simplest lattice Boltzmann equation (LBE) is the lattice Bhatnagar-Gross-Krook (BGK) equation, 
based on a single-relaxation-time approximation. Due to its extreme simplicity, the lattice BGK (LBGK) 
equation has become the most popular lattice Boltzmann model in spite of its well known deficiencies, such 
as the instability at low viscosity. 

We have constructed two three-dimensional LBE models with 15 and 19 velocities based on the multiple- 
relaxation-time (MRT) lattice Boltzmann equation. The MRT LBE model (also referred to as the generalized 
lattice Boltzmann equation (GLBE) or the moment method) overcomes some obvious defects of the LBGK 
model, such as fixed Prandtl number (Pr = 1 for the BGK model) and fixed ratio between the kinematic and 
bulk viscosities. The MRT lattice Boltzmann equation has been persistently pursued, and much progress has 


19 


been made in the following areas: formulation of optimal boundary conditions, interface conditions in multi- 
phase flows and free surfaces, thermal and viscoelastic models, models with reduced lattice symmetries, 
and improvement of numerical stability. It should be stressed that most of the above results cannot be 
obtained with the LBGK models. Applying optimization techniques in coding, the computational efficiency 
of the RLBE method can be fairly close to that of the LBGK method for most practical applications (RLBE 
schemes could be ca. 15% slower than their LBGK counterparts in terms of the number of sites updated 
per second). Recently it was shown that the multiple-relaxation-time LBE models are much more stable 
than their LBGK counterparts, because the different relaxation times can be individually tuned to achieve 
‘optimal’ stability. 

The present work is a result of collaboration with D. d Humieres (Ecole Normale Superieure, France). 
I. Ginzburg (Fraunhofer-Institut fiir Techno- und Wirtschaftsmathematik, Germany), M. Krafczyk (Teehnis- 
che Universitat Braunschweig, Germany), and P. Lallemand (Laboratoire ASCI-CNRS, Universite Paris-Sud 
(Paris XI Orsay), France). The present work has been partially funded by the United States Air Force Office 
for Scientific Research under Grant No. F49620-01-1-0142 (technical monitor: J. Tishkoff). 

We intend to apply the LBE MRT model for DNS and LES simulations. 
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STRUCTURES AND MATERIALS 


Wholly Aromatic Ether-imide Liquid Crystals I. Synthesis and Characterization 

Thao Ding emam 

In order to improve the processability of single- wall and multi-wall carbon nanotubes, we have synthesized 
a series of wholly aromatic ether-imide model compounds composed of an electron-poor imide core and 
electron-rich aryl-ether flexible tails. We choose to work with this class of compounds because aromatic imides 
are able to form charge transfer complexes with closed-shell carbon nanostructures, and this interaction can 
be utilized to solubilize carbon nanostructures. 

To date we have synthesized and characterized a novel series of low-molar-mass ether-imide rod-shaped 
model compounds. All ether-imides were obtained by terminating the appropriate rigid core dianhydride, i.e., 
pyromellitic dianhvdride (PMDA), 1,4,5,8-napthalenetetracarboxylic dianhydride (NDA), 3,3\4,4’-biphenvl- 
tetracar boxy lie dianhydride (BPDA), and 3,3’,4,4’-oxydiphthalic dianhvdride (ODPA) with three different 
flexible aryl-ether tails. The mono-functional aryl-ether amines, i.e., 4-(3-phenoxy-phenoxy)-phenvlamine 
and 4-(3-phenoxy-3-phenoxy-phenoxy)-phenylamine, were synthesized using standard fluoro-displacement 
and Ullmann condensation techniques. The corresponding ether-imide model compounds were obtained 
in high yields using a one-step solution imidization procedure. Increasing the number of meta - substituted 
aryl-ether units reduces the melt transition temperatures and at the same time it increases the solubility of 
the model compounds. Most model compounds are crystalline solids and form isotropic melts upon heating. 
2,7-Bis-(4-phenoxy-phenyl)benzo[/mn][3,8]phenanthroline-l,3,6,8-tetraone, however, displays a smectic A and 

rectangular columnar phase when cooled from the isotropic phase. This is the first example to date, in which 
a mesophase is detected in a wholly aromatic ether-imide compound. Increasing the number of flexible me.ta- 
substituted arvl-ethers results in low melting ether-imides that super-cool and form amorphous glasses with 
well defined glass-transition temperatures. 

Future work will include temperature-dependent NMR, Raman spectroscopy, and atomic force mi- 
croscopy (AFM) experiments, which will allow us to monitor and quantify the interactions between the 
model compounds and carbon nanotubes. In addition, we will use molecular modeling as a tool to predict 
which imide model compounds are the most promising candidates for forming charge transfer complexes 
with carbon nanotubes. 

This research was conducted in collaboration wit h Kris Wise and Cheol Park (ICASE), and Peter Lillehei 
(NASA Langley). 

Carbon Nanotube Pullout from a Polymer Matrix in a Nanocomposite 

Sarah-Jane V. Frankland and Vasyl Michael Harik 

Carbon nanotubes have been proposed as fillers to reinforce polymer matrices. The ability to transfer load 
across the polymer-nanotube interface is fundamental to the mechanical performance of these materials. To 
characterize the nanotube-polymer interfacial properties, a detailed theoretical study of the nanotube pullout 
from polymer matrices is being carried out. The deformation processes occurring in the composite as the 
nanotube moves within the matrix are relevant to the design and fabrication of polymer-nanotube composites 
as w r ell as interfacial failure mechanisms. The objective of this study is to use atomic-level information from 
molecular dynamics simulation to identify deformation processes during pullout, to parameterize the related 
shearing behavior from the simulations, and to develop higher-level models that incorporate the simulated 
input. 
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Molecular dynamics simulations have been carried out for the process of a nanotube pullout from a 
crystalline polyethylene matrix. Semi-periodic variations in the nanotube velocity are observed which can 
be related to the local structure of a nanotube. Analysis of the force- velocity dependence has led to being 
able to fit the simulation data into a model. This model gives an estimate of the local effective viscosity. In 
this way, the shearing deformation of a nanotube-polymer nanocomposite may be evaluated. 

Simulation data will be generated for non-bonded and chemically functionalized polymer-nanotube com- 
posites to further identify and model processes during nanotube pullout. Analogous viscosity information 
will be generated for and applied to the slippage of carbon nanotubes within a bundle in the design of woven 
nanotube fibers. 

This work was done in collaboration with J.A. Hinkley and T.S. Gates (NASA Langley). 

Transverse Moduli of Carbon Nanotube Ropes 

Sarah-Jane V. Frankland 

Single-walled carbon nanotubes occur in hexagonal packed bundles or ropes. The mechanical properties 
of these bundles as groups of nanotubes rather than individual nanotubes is not well-characterized empirically 
or theoretically. In particular, the transverse properties which depend on the weak non-bonding van der 
Waals interactions between the nanotubes have been generally ignored in favor of the more impressive 
longitudinal properties of perfect carbon nanotubes such as the Young’s modulus typically reported to be 1 
TPa. The transverse properties must be considered in the design of polymer composites or other materials 
such as woven nanotube yarns, which depend on nanotubes for reinforcement. The objective here is therefore 
to assist the design effort of these materials by providing reasonable computational estimates of the transverse 
moduli of carbon nanotube bundles. 

In this work, the transverse moduli of carbon nanotube bundles have been calculated from a simple 
molecular statics model of rigid hexagonally packed nanotubes about 1 nm each in diameter. The intertube 
potential depends solely upon the van der Waals interactions between the nanotubes which are represented by 
the Lennard-Jones potential. The rigid approximation is reasonable for nanotubes of this size. Calculations 
have been done for the elastic constants C22, E22, bulk modulus K23, Poisson’s ratio v23, and shear modulus 
G23, where 2 and 3 are the directions transverse to the nanotube axis. These independently determined values 
relate consistently via the Hash in relations. Comparisons were made with existing data in the literature. 
These calculations provide an estimate for the transverse behavior of nanotube bundles based on relatively 
simple consideration of van der Waals interactions. 

This work will be extended to the calculation of other mechanical properties for additional nanotube 
ropes. Refinement of the nanotube-nanotube interaction potential is being considered. The moduli will be 
used in designing and evaluating the mechanical properties of nanotube yarns and polymer composites. 

This work was done in collaboration with J.A. Hinkley, T.S. Gates, and E. Saether (NASA Langley) 
and R.B. Pipes (University of Akron). 

Mechanics of Carbon Nanotubes 

Vasyl Michael Harik 

Carbon nanotubes possess extraordinary physical properties (e.g., high stiffness and strength along with 
enormous electrical and thermal conductivities). Potential applications range from new electronic devices and 
nanotube-based scanning probes to multifunctional polymer films and sensors for new aerospace structures. 
The goal of this research is to develop constitutive models for carbon nanotubes and a methodology for the 
implementation of continuum models at the nanometer scale. 
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Ranges of validity and length-scale limitations of the continuum models for carbon nanotubes under 
compressive loads arc* examined. Hierarchical dimensional analysis of the nanomechanical buckling problem 
and the geometric- and material parameters involved is carried out. As a result, three c lasses of carbon 
nanotubes are identified. One of the classes consists of the newly predicted carbon nanobeams. The scaling 
laws and the key non-dimensional parameters controlling the deformation process of nanotubes have been 
derived. Expressions for critical strains are tailored for various nanotubes and their values are compared 
with the molecular dynamics simulations. Model applicability maps are constructed for different nanotube 
geometries in the parameter space. 

Future research involves analysis of the nanotube-based AFM tips and the effects of nanotube geometry 
on the; constitutive properties of nano-structured materials. Applicability of advanced shell theories will be 
examined as well. 

This research is conducted in collaboration with T.S. Gates, M.P. Nemeth, and D.R. Ambur (MDB. 
NASA Langley). 

Modeling of Multifunctional Composite Materials 

Vasyl Michael Harik 

Design of new aerospace structures, such as flexible wings, structures with active controls, multifunc- 
tional membranes, and inflatable antennas require fundamental understanding of electromechanical coupling 
in material behavior so that they can be modeled accurately in structural analyses. Nanostructured mate- 
rials offer additional benefits that stem from their superior physical properties. The goal of this research 
is to develop hierarchical constitutive models that link the nano-scale molecular structures with continuum 
micromechanics. Such models are critical for the development, design, and optimization of novel nanocom- 
posites. 

Constitutive models for multiphysics phenomena are examined along with the structure of carbon- 
filled polymer matrix systems. Homogenization criteria for hierarchical connectivity and thermodynamically 
consistent material averaging are derived for nanotube-reinforced polymer composite systems. This modeling 
approach has been used to link molecular dynamics simulations of nanocomposites with nanomechanical 
models that are analogous to similar models in micromechanics. 

The future work is directed toward a better understanding of applicability and limitations of this method- 
ology and alternative techniques for multi-scale modeling. 

This research is carried out in collaboration with S.J.V. Frankland (ICASE), D.R. Ambur and T.S. 
Gates (MDB, NASA Langley). 

Stitched Stiffened Panels Loaded in Shear 

Navin Jannky 

Progressive failure studies of composite panels were carried out for stiffened stitched panels. These studies 
were carried out to validate modeling methodologies for intralaminar damage modes being developed. The 
development of such damage modeling methodologies are discussed in earlier reports. 

Intralaminar damage modes such as matrix cracking, fiber-matrix shear, and fiber failure are modeled 
by degrading the material properties. The stitched panels were stiffened panels and were loaded in in-plane 
shear. Two stiffened panels were considered, one had a notched and the other was un-notched. These panels 
were loaded well into their postbuckling regime. Experimental results from the test include strain field data 
from video image correlation in three dimensions (VIC-3D) in addition to other strain and displacement 
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measurements. Progressive failure analyses of these panels did not include interlaminar damage mode since 
the panels are stitched. Initial geometric imperfections and thickness imperfections were included in the 
finite element models. The three intralaminar failure modes considered for the stitched panels in the study 
accurately represent the damage scenario in the postbuckling regime. The analytically determined panel 
response, failure modes, and damage locations compare well with the experimental results for the pointwise 
as well as for full field measurements. Including thickness imperfection provides better agreement with the 
experimental results. Progressive failure analyses results with successive loading and unloading sequences 
to mimic the experimental loading conditions are also in good agreement with the experimental results. It 
was found that these progressive failure analyses can account for the damage accumulated in a previous run. 
Damage detection after failure was attempted on the un-notched panel using a thermographic scan. 

Future work will focus on more comparison of progressive failure analyses for different types of loading. 
Damage detection at different load levels during the experiments will be attempted using thermographic 
scan. Full field measurement using VIC-3D will also be used. Attempt to incorporate physics based failure* 
criteria in the progressive failure analyses will be made. 

This work is done in collaboration with Damodar R. Ambur and Mark W. Hilburger (NASA Langley). 

Intralaminar and Interlaminar Progressive Failure Analyses 

Navin Jaunky 

A progressive failure methodology is developed that includes the interaction of intralaminar and inter- 
laminar damage modes. The intralaminar failure is as described above. The interlaminar failure mechanism 
such as delamination is simulated by positioning interface elements between adjacent sublaminates. A non- 
linear constitutive law is postulated for the interface element that accounts for a multi-axial stress criteria 
to detect the initiation of delamination, a mixed-mode fracture criteria for delamination progression, and 
a damage parameter to prevent restoration of a previous cohesive state. The methodology is validated us- 
ing experimental data available in the literature on the response and failure of quasi-isotropic panels with 
centrally located circular cutouts loaded into the postbuckling regime. 

Very good agreement between the progressive failure analyses and the experimental results is achieved if 
the failure analyses includes the interaction of intralaminar and interlaminar failures. The analyses demon- 
strated the following: 

1. The intralaminar damage induced the initiation of delamination at the free edge of the cutout at 
the same locations identified in the experiments. 

2. Delamination growth occurred perpendicular to the loading direction in regions where large com- 
pressive loads were present. 

3. Catastrophic failure was a consequence rapid delamination growth. 

4. Delamination growth followed by delamination arrest leads to a residual strength, a finding consistent 
with the experimental data. After further loading, delamination growth occurs leading leading to 
catastrophic failure. 

Future work will focus on more comparisons of progressive failure analyses for different types of loading. 
Damage detection at different load levels during the experiments will be attempted using thermographic 
scan. Full field measurement using VIC-3D will also be used. Attempt to incorporate physics based failure 
criteria in the progressive failure analyses will be made. 

This work is done in collaboration with Vinay Goyal and Eric Johnson (Virginia Polytechnic Institute 
and State University). Damodar R, Ambur (NASA Langley) was also involved. 
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Debonding in Stringer Reinforced Composite Components 

Ronald Krueger 

Many composite components in aerospace structures are made of flat or curved panels with eo-cured or 
adhesively bonded frames and stiffeners. Previous investigations of the failure of secondary-bonded st ruct ures 
focused on loading conditions typically experienced by aircraft crown fuselage panels. First tests were 
conducted with specimens cut. from a full-size panel to verify the integrity of the bondline between the 
skin and the flange or frame. A simpler and cheaper specimen configuration that would allow detailed 
observations of the failure mechanism at the skin/flange interface was proposed later. The investigations 
focused on the failure mechanisms of a bonded skin/flange coupon configuration subjected to tension, three 
and four-point, bending, and combined tension/bending loading. The failure that initiates at the tip of the 
flange in these specimens is identical to the failure observed in the full-scale panels and the frame; pull-off 
specimens. An analytical methodology was developed to predict the location and orientation erf the first 
transverse matrix crack based on the principal transverse tension stress distribution in the off axis plies 
nearest the bondline in the vicinity of the flange tip. A fracture mechanics approach was used to investigate 
delamination onset once the initial crack had formed. The initial crack was modeled as a discrete discontinuity 
at. a location suggested by the microscopic investigation. The complex nature of the failure observed during 
the experiments, where the delamination changed across the specimen width from a delamination running 
at the skin surface 45°/-45° layer interface to a delamination propagating in the bondline above, suggests 
the need for a three-dimensional model. Since many layers of brick elements through the thickness would 
be required to model the individual plies, the size of three-dimensional finite element models, however, may 
become prohibitively large. Two-dimensional models of a longitudinal cut through the specimen allow for a 
detailed modeling of the individual plies and the adhesive in thickness direction. Two-dimensional models are 
preferred by industry due to the fact that modeling time, as well as computational time, remains affordable, 
especially if many different configurations have to be analyzed during the initial design phase. However, 
it is inherent to any two-dimensional finite element model that the geometry, boundary conditions, and 
other properties are constant across the entire width. The effect of two-dimensional modeling assumptions 
is most marked for 45° plies because of their high in-plane Poisson’s ratio, while it is small for 0° and 90°. 
The objective of the current, study was to evaluate how the results obtained from two-dimensional finite 
element models of the specimen compared to data obtained from full three-dimensional simulat ions in order 
to find the limitations of using simple two-dimensional models. Investigating the limits and usefulness of 
any numerical model is in general an important engineering science problem that reaches beyond the current 
analysis of the skin / stringer debond specimen. 

Finite Element analyses were performed for simple skin/flange specimens subjected to tension and bend- 
ing loads. Geometrically nonlinear finite element analyses using two-dimensional plane-stress and plane- 
strain elements as well as three different generalized plane strain type approaches were performed. The 
computed skin and flange strains, transverse tensile stresses and energy release rates were compared to 
results obtained from three-dimensional simulations. The plane stress and plane strain models provided 
results for skin and flange strains, as well as energy release rates, which form an upper and low r er bound of 
the results obtained from full three-dimensional analysis. The results from generalized plane strain models 
fall between the results obtained from plane stress and plane strain models. The generalized plane strain 
models capture the surface strains on the flange and skin very well. Computed energy release rates art; 
within 9% of the results from three-dimensional analysis, however, there is not consistency across the load 
cases with respect to which model performs best. Stresses obtained from plane strain and generalized plane 
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strain models, were excessively high when compared to results from three-dimensional analysis. The stresses 
from full three-dimensional analysis were lower than results from any of the* two-dimensional models, but arc' 
closest, to the plane stress results. Based on the results of this investigation it is recommended to use results 
from plane stress and plane strain models as upper and lower bounds. Two-dimensional models may also bo 
used to qualitatively evaluate the stress distribution in a ply and the variation of energy release rates and 
mixed mode ratios with delamination length. The current recommendations are based on the analysis of the 
skin /stringer specimen with a particular layup. Further analyses are required before the recommendations 
may be generalized. 

The current study presents an intermediate step where three-dimensional models were created by extrud- 
ing two-dimensional models across the width. The fact that the delamination changed across the specimen 
width from a delamination running at the skin surface 45°/-45° layer interface to a delamination propagating 
in the bondline above, however, is still not accounted for in this model. Nevertheless, the three-dimensional 
model takes width effects into account and therefore provides insight into the limitations of the use of two- 
dimensional finite element models. For future detailed modeling and analysis of the damage observed during 
the experiments, the shell/3D modeling technique offers great potential for saving modeling and computa- 
tional effort because only a relatively small section in the vicinity of the delamination front needs to be 
modeled with solid elements. 

This work is done in collaboration with Isabelle L. Paris (National Research Council), T. Kevin O’Brien 
(Army Research Laboratory, Vehicle Technology Directorate located at NASA Langley), and Pierre J. 
Minguet (The Boeing Company, Philadelphia). 

Intelligent Optics 

Mark Little 

Intelligent optics, materials, or devices whose optical properties may be controlled by external means, 
have a significant role in many areas of today’s society. Some include the communications industry, scientific 
research, and space applications. One problem with current intelligent optical materials is a characteristically 
slow response time. For use in communications and space it is imperative to have fast optical devices. One 
possible solution is to exploit the Stark and Zeeman effects in materials. These quantum mechanical lifting 
of electronic degeneracy by electric or magnetic fields may induce a change in the materials optical properties 
that can be exploited. Since these transitions are quantum in nature their response times will be very fast. In 
addition, current optics are designed to work at a set frequency range. If materials could be developed that 
would allow the frequency range of standard optical devices to be changed without the need to continually 
change components, the optics world would be revolutionized. The goal of this research is to fabricate new 
optical thin film (OTF) materials that use the Stark, Zeeman, or other effects as its active optical ‘switch’. 

Current work is focused on the exploration growth for new OTF materials. The growth systems for the 
lab have been installed and are producing products. The Stark and Zeeman effect will require an atomically 
isolated dopant level residing in a host material. Many dopant candidates are being tried and the optical 
properties with applied fields are being measured. Various host materials, including standard optical glasses, 
oxides, wide band gap semiconductors (both crystalline and amorphous), and ferroelectric materials are 
being investigated. In addition to the optical characterization, basic electronic testing capabilities of these 
materials have been established to check their suitability as new electronic materials. 

Future work is to include the continued growth and characterization of new materials for smart optical 
devices based on the Stark and Zeeman effects. Also new methods for creating active diffraction gratings are 
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in the infancy stage. Brainstorming is occurring to determine a feasible 1 method of achieving this goal. Then 1 
is interest in this area of research from the Army Office of Research and a collaborative effort is attempting 
to be established. I would like to thank Glen King and Sang Choi for their help and guidance with this 
research. 

Constitutive Modeling of Nanotube/Polymer Composites 

Gregory Odegard 

In the last few years, SWNT-reinforced polymer composites have generated considerable interest in the 
materials research community partly because of the potential for large increases in strength and stiffness 
when compared to typical carbon-fiber-reinforced polymer composites. In order to facilitate the development 
of nanotube-reinforced polymer composites, constitutive relationships must be developed to predict the bulk 
mechanical properties of the composite as a function of molecular structure and bonding of the polymer, the 
nanotube, and the polymer /nanotube interface. 

A modeling methodology has been developed that predicts continuum mechanics-based elastic constants 
of nanotube composites using atomic-scale information. Molecular structure and bonding information, which 
is obtained through Molecular Dynamics simulations, is used to construct an equivalent-continuum model of 
a nanometer- scaled representative volume element of the nanotube composite. Specifically, models have been 
generated for nanotube/polyimide and nanotube/ Polypropylene composites (both with a PmP\ interface). 
It has been shown that the predicted elastic properties agree with other multi-scale modeling methodologies 
and with experimental results. Therefore, the model is ready to be implemented as a tool in the development 
of nanotube/polymer composites. 

Future work will involve the modeling of additional nanotube composite 1 systems. In order to under- 
stand the effect that covalent bonding between the nanotube and polymer has on the overall mechanical 
properties of the composite, functionalized and non-functionalized nanotube/polyethylene composites will 
be modeled. Nanotube/liquid crystal polymer systems will be modeled to explore their potential use as 
structural materials. 

This research has been performed in collaboration with T.S. Gates, E.J. Siochi (NASA Langley), S.J.\ . 
Frankland, V.M. Harik, C. Park, K.E. Wise (ICASE), R.B. Pipes (University of Akron), and P. Hubert (Old 
Dominion University). 

Nanoindentation of Nano-structured Materials 

Gregory Odegard 

Nanoindentation is a promising method of measuring mechanical properties of materials at the nanometer- 
length scale. Among the advantages of nanoindentation is the ability to examine surface effects and individual 
constituents of a structural material. Because nanoindentation has been used almost exclusively for metallic 
and ceramic materials, its application to characterizing viscoelastic properties of nano-structured polymers 
and polymer composites needs to be explored in detail. 

A series of tests has been performed on various advanced structural polymers for measurement of vis- 
coelastic properties. The dependence of measured properties on testing parameters (e.g., strain rate, har- 
monic displacement, and harmonic frequency) has been investigated. Procedures for mounting bulk samples 
and thin films have been developed. Comparisons of measured viscoelastic properties obtained from other 
test methods show remarkable agreement. Together, these studies have established an overall testing proce- 
dure that can be used for testing a wide range of polymers under various test conditions. 
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Future work will involve the viscoelastic characterization of nanotube/polymer composite films and other 
nano-structured polymeric materials. 

This research has been performed in collaboration with T.S. Gates (NASA Langley). 

Corona Poling of Partially Cured Polyimide 

Cheol Park 

In-situ poling and imidization of the partially cured (/?-CN)APB/ODPA was studied in an attempt to 
maximize the degree of dipolar orientation and the resultant piezoelectric response. The dielectric relaxation 
strength, remnant polarization, and piezoelectric responses were evaluated as a function of the poling profile. 
The partially cured, corona-poled polymers exhibited higher dielectric relaxation strength (<5c), remnant 
polarization (P r ), and piezoelectric strain coefficient (<i 3 3) than the fully cured, conventionally poled ones. 

Application of the in-situ poling and imidization of partially cured polyimides for amorphous polvimides 
containing higher dipolar concentration will be investigated. 

This research was performed in collaboration with Zoubeida Ounaies (Virginia Commonwealth Univer- 
sity), Kristopher E. Wise (ICASE), and Joycelyn S. Harrison (NASA Langley). 

Single Wall Carbon Nanotube Polymer Composites 

Cheol Park 

Carbon nanotubes are of great interest because of their unique electronic and mechanical properties 
combined with their chemical stability. Single wall carbon nanotubes (SWNTs), however, have been rarely 
used as an electrical or mechanical inclusion in a polymer matrix mainly because of the difficulty of efficient 
dispersion. A novel process to effectively disperse SWNTs in an aromatic polymer was developed. This 
process involves an in-situ polymerization of monomers of interest in the presence of sonication during the 
polymerization process. The goal of this study is to develop a method to disperse SWNTs into polymer ma- 
trices on a nanoscale level. The resultant SWNT-polymer nanocomposite exhibited significant conductivity 
enhancement (10 10 order of magnitude) at a very low loading (0.1 wt%) without sacrificing significant optical 
transmission. Mechanical properties as well as thermal stability mechanical properties were also improved 
with the increase of the SWNT incorporation. The SWNT-polymer nanocomposites are useful in aerospace 
and terrestrial applications. 

SWNT/polymer nanocomposites will be further characterized microscopically (SPM, HRTEM, and 
HRSEM). A combined experimental and theoretical study will be performed to understand the nature of 
the interaction and identifying new' complexes with improved properties. 

This research w^as performed in collaboration with Zoubeida Ounaies (Virginia Common wealth Uni- 
versity), Sarah J. Frankland (ICASE), Kristopher E. Wise (ICASE), Greg Odegard (ICASE), Torn Gates 
(NASA Langley), Sharon E. Low'ther (NASA Langley), Peter Lillehei (NASA Langley), Mia Siochi (NASA 
Langley), Joycelyn S. Harrison (NASA Langley), Alice Chang (Lockheed- Martin), Roy Crooks (Lockheed- 
Martin), Neal D. Evans and Edward A. Kenik (Oak Ridge National Laboratory), Jim Bentley (ORNL). and 
Satish Kumar (Georgia Tech). 

Scanning Electron Microscopy of Carbon Nanotube Materials 

Jason Rouse 

A central difficulty in analyzing the quality of polymer/carbon nanotube composites is the ability to 
accurately determine the size of the nanotube materials in the composite, the uniformity of their dispersion. 
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and their relative orientation. While relative average orientation can be probed by polarized Raman spec- 
troscopy and optical microscopy can be used to identify very poorly dispersed nanotube samples, the ability 
to characterize highly dispersed nanotube samples has mostly relied on transmission electron microscopy 
(TEM). A major drawback of TEM is the time consuming process of preparing highly uniform ultrathin (15 
nm ) samples for imaging. With the advent of new high-resolution scanning electron microscopes, the ability 
to image nanotubes in the polymer matrix without time consuming sample preparation is being investigated. 

Currently, polymer/carbon nanotube composites with varying nanotube loading levels, different polymer 
matrixes, and produced under multiple processing conditions are being imaged. Early results indicate that 
small changes in imaging parameters can dramatically affect the quality of the images obtained. 

This research is performed in collaboration with P.T. Lillehei and Denis Working (NASA Langley), and 
Cheol Park (ICASE). 

Stepwise Assembly of Polymer/Carbon Nanotube Multilayered Films 

Jason Rouse 

An area of difficulty in forming polymer/carbon nanotube composites is the ability to sufficiently disperse 
the carbon nanotubes in the polymer matrix. Additionally, at high nanotube loading the high viscosity of the 
resulting polymer/carbon nanotube melt makes processing difficult. In order to overcome these difficulties, 
we have focused on the stepwise assembly of polymer/carbon nanotube films using molecular interactions. By 
forming films at the nanometer level, it is possible to investigate film structures that can not be constructed 
using bulk processing techniques. 

As a starting point for the formation of polymer/single- walled carbon nanotube films, we explored 
the use of electrostatic interactions. The formation of thin, uniform films by the sequential adsorption of 
polyionic species of opposite charge is an area of significant research interest. We have found that the 
oxidative procedure used to remove amorphous carbon during the purification of carbon nanotubes also 
results in the formation of oxidized defects on the tubes themselves. These oxidized groups allow carbon 
nanotube bundles to act as polyionic species. Atomic force and scanning electron microscopies revealed that 
the sequential absorption of the cationic poly(diallyldimethylammonium chloride) followed by single walled 
carbon nanotubes results in the formation of dense films of overlapping nanotubes. 

Future work will involve the characterization of the mechanical and electrical properties of these films. 
Additionally, replacement of poly(diallyldimethylammonium chloride) with a more thermally and mechani- 
cally' stable polymer will be explored. 

This work was done in collaboration with P.T. Lillehei (NASA Langley). 

Formation of Novel Organic-inorganic Electroactive Composites 

Jason Rouse 

The rigors of low earth orbit require materials that possess a unique combination of physical prop- 
erties. Resistance to atomic oxygen and ultraviolet radiation and the ability to tolerate large temper- 
ature changes are required. A material that possesses these attributes is the inorganic-organic polymer 
poly(dimethlysiloxane) (PDMS). Unfortunately, PDMS lacks functional groups that, respond to electrical 
fields. The incorporation of such groups into PDMS is the focus of this research. 

Novel PDMS-based hybrid composites were prepared by the tin-catalyzed condensation of cyanoethvl- 
triethoxysilano with hydroxy-terminated PDMS oligomers. Films prepared at low catalyst concentration 
resulted in the phase-separation of the silicate into micron-sized spherical particles. These films were flexi- 
ble, thermal robust, and responded to electrical fields. Measurements revealed that the dielectric constant 
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appeared to be frequency independent in the range lOOHz-lOOkHz, possibly the result of the phase separa- 
tion. Electrostatic forces dominated the strain mechanism, where the overall strain response was proportional 
to E2, and was dictated by the dielectric constant of the composite and the driving voltage. 

We are currently preparing composites with a higher loading of catalyst to retard phase separation. 

This research is conducted in collaboration with Zoubeida Ounaies (Virginia Commonwealth University). 
Mia Siochi (NASA Langley), and Peter Lillehei (NASA Langley). 

Low-color Polyimides for Space Applications 

Kent A. Watson 

NASA has several space applications that currently use or are evaluating polyimides. These include thin 
films as membranes on antennas, concentrators, coatings on second-surface mirrors, solar sails, sunshades, 
thermal/optical coatings and multi-layer thermal insulation blanket materials. Depending on the application, 
the film will require a unique combination of properties. These may include atomic oxygen resistance, 
ultraviolet and vacuum ultraviolet radiation resistance, low color /low solar absorption and high mechanical 
properties (strength, modulus, and toughness). 2,3,3\4’-biphenyltetracarboxylic dianhvdride (a-BPDA) was 
investigated as a dianhydride for the synthesis and characterization of a series of polyimides. Polyimides 
derived from a-BPDA were thought to offer an attractive combination of properties such as low color, 
good solubility, high thermal emissivitv, low solar absorptivity, and high tensile properties. This work was 
performed to obtain fundamental information on chemical structure/property relationships that could lead 
to the development of improved films for space use. 

We successfully prepared a series of polyimides derived from a-BDPA and compared the materials with 
the corresponding 3,3’4,4’-biphenyltetracarboxylic dianhydride (s-BPDA). a-BPDA was less reactive with 
aromatic diamines than the s-BPDA. Polyimides from a-BPDA had higher glass transition temperatures, 
less color (higher optical transparency), and lower thin film tensile properties than polyimides derived from 
s-BPDA. A polvimide derived from a-BPDA and l,3-bis(3-aminophenoxy)benzene (APB) had the lowest 
color of all polyimides prepared. Advantage of this material when compared to state-of-the-art polyimides 
currently in use is the potentially much lower cost of the material. 

Future work includes continued evaluation of the polyimide prepared from a-PBDA and APB for po- 
tential use as a low-color material for space applications. Future work also includes the incorporation of 
single-wall carbon nanotubes into the most promising materials in order to provide sufficient conductivity 
to mitigate static charge build-up. 

This work was performed in collaboration with Paul Hergenrother, Joseph Smith. Jr., and John Connell 
(NASA Langley) and R. Yokota (Institute of Space and Astronautical Science). 

Organic Nanocomposites from Functionalized Carbon Nanotubes 

Kent A . Watson 

Thin film polymer materials represent enabling technology for future spacecraft. These spacecraft include 
large ultra-lightweight space structures, space-based observatories, and solar sails. Currently no materials 
exist which possess the desired combination of properties necessary to perform in these missions. Spacecraft 
in low Earth orbit (LEO) typically have lifetimes of five to ten years and during this time the materials 
are bombarded with ultraviolet (UV) arid vacuum ultraviolet (VUY) radiation, atomic oxygen (AO), mi- 
crometeroids, proton radiation, and electron radiation. Materials have recently emerged that show improved 
resistance to AO, UV, and VUV radiation. However, all of these new materials are insulators (i.e., non- 
conducting). Due to charged species being present in the space environment, these polymers can build up 
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significant electrical charges that can eventually lead to sparking or electrical arcing. Transfer of the electri- 
cal charge from the film-based material to the main body of the spacecraft can cause catastrophic* damage 
to computers and sensitive electronic devices. Recent work has involved incorporation of single- wall carbon 
nanotubes (SCNT) into polymeric materials to add sufficient conductivity to mitigate static charge' build-up. 
One main obstacle to using SCNTs is their difficulty in dispersing in an organic medium. They are both 
insoluble in organic solvents and the SCNTs tend to agglomerate into bundles. To overcome these difficulties 
we will investigate the preparation of organic nanocomposites using functionalized SCNTs. 

We arc* currently investigating different methods of functionalizing SCNTs primarily by taking advantage 
of functional groups already existing on the surface of the nanotubes. Organic functional groups are formed 
on the surface of the nanotubes during purification (acid treatment). Oligomer or polymer chains can be 
chemically bound to the surface of the nanotubes at these functional group sites. Attachment of oligomer or 
polymeric organic moieties should aid in dispersion of the SCNTs in organic solvents and also show improved 
compatibility with host polymers. Our preliminary work has involved the preparation of functionalized 
SCNTs and incorporation into space durable polymers. We have prepared functionalized SCNTs that contain 
small molecules as well as oligomeric groups attached to the nanotubes surface. 

Current and future work includes exploring additional methods for SCNT functionalization and char- 
acterization of the functionalized SCNTs. Characterization of the polymer/functionalized SCNT nanocom- 
posites will also be investigated. 

This work is being performed with the collaboration of John Connell (NASA Langley) and Sun 
(Clemson University). 

Modeling Nanotube Dispersion in a Polymer Matrix 

Kristopher Wise 

A key enabling technology for the application of carbon nanotubes in practical devices is the ability to 
disperse them in a polymer matrix. Currently, kinetically stable (short lived) nanotube-polymer solutions 
can be prepared by a variety of traditional stirring/mixing methods. Unfortunately these solutions eventually 
phase separately due to the strength of the tube-tube (relative to the polymer-tube) interactions. Achiev- 
ing thermodynamically stable nanotube-polymer composite systems will require enhancing the attractive 
interaction between the matrix and nanotubes. 

Using a combination of high-level ab initio (Hartree-Fock and Density Functional Theory) and classical 
(Molecular Dynamics and Monte Carlo) simulation methods, we are studying the effect of varying the 
chemical composition of the polymer on its ability to stabilize nanotube dispersions. Based on preliminary 
results, we have proposed a mechanism involving the formation of an electron donor/acceptor complex 
between particular functional groups on the polymer and the nanotube. 

Ongoing work is proceeding in two primary directions: expanding the types of polymer systems being 
studied and developing new simulation methodologies. Based on what has been learned to this point, we 
have a number of new polymer systems which we will be testing both experimentally and through simulation. 
We are also working on ways of expanding our simulation capabilities to encompass larger composite systems 
and longer simulation times. 

This work is being performed in collaboration with Cheol Park (ICASE), Zoubeida Ounaies (Virginia 
Commonwealth University), and Mia Siochi and Joycelyn S. Harrison (NASA Langley). 
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Development of Facilities for Characterizations of Electromechanical Properties of 
Electroactive Polymer-based Materials and Devices 

Tian-Bing Xu 

Electroactive polymers (EAP), which change shape as an electric field is applied, can he used in many 
areas such as artificial muscles and organs, smart materials for vibration and noise control, electromechanical 
actuators and sensors for robots, acoustic transducers used for underwater navigation and medical imaging, 
and fluid pumps and valves. Compared with piezocerarnic and magnetostrictive actuator materials, polymers 
have many advantages such as flexibility, easy processing, light weight, arid low cost. Polymeric materials 
can also withstand a large dimensional change (strain) without fatigue and are quite robust. Recently, a 
Graft elastomer-based EAP with giant strain, relatively high elastic modulus, and high energy density has 
been developed at NASA Langley. There are, however, no commercialized methods for electromechanical 
property characterizations in the EAP films. In order to keep world leadership in the field of research on 
EAP and EAP-based devices, it is necessary to develop/install the unique world leader facilities at NASA 
Langley to match the needs of characterization of electromechanical response on the EAP films and EAP 
based electromechanical devices. 

A transverse strain measurement set-up (Strain response in lateral direction of the films), a longitudinal 
strain set-up (Strain response in the thickness direction), and a force measurement set-up for electrome- 
chanical response in EAP films have been designed and fabricated in this duration. The performances of 
these facilities are under evaluation. On the other hand, a four probes conductivity measurement station has 
been installed and evaluated. It works fine. In addition, the laser vibrometer based displacement/actuation 
measurement set for micromachined actuators has been designed and the laser vibrometer is selected (un- 
der order process). There unique facilities can greatly increase our capability to do/lead research on EAP 
materials and EAP-based devices. 

We are going to characterize the electromechanical response of EAP films at different process conditions 
and compositions to optimize/improve the electromechanical properties of the EAP films after calibrating 
the new facilities. On the other hand, we also plan to develop the EAP-based actuators and micro-electro- 
mechanical systems (MEMS) for aerospace applications. 

This work is accomplished in collaboration with Ji Su and Joycelyn S. Harrison (NASA Langley). 


32 


I 



REPORTS AND ABSTRACTS 


Stals, Linda: The solution of radiation transpor t equations with adaptive finite elements. ICASE Report No. 
2001-29 , ( N A S A / C R - 2 00 1 -2 1 1 230 ) , October 15. 2001, 20 pages. Proceedings of the Tenth Copper Mountain 
Conference on Multigrid Methods. 

We compare the performance of an inexact Newton-multigrid method and Full Approximation Storage 
rnultigrid when solving radiation transport equations. We also present an adaptive refinement algorithm and 
explore its impact, on the solution of such equations. 

Shin, Jong-Yeob: Analysis of linear parameter varying system models based on reachable sets. ICASE Report 
No. 2001-30 , (NASA/CR-2001-211231), October 24, 2001. 22 pages. To be submitted to American Control 
Conference (ACC) 2002. 

This paper presents the analysis method of quasi-LPV models, comparing the ellipsoid set which contains 
the reachable set of a nonlinear system to define which quasi-LPV model is less conservative to represent the 
nonlinear dynamics. Three quasi-LPV models are constructed from a nonlinear model using three different 
methods, to facilitate synthesis of an LPY controller for the nonlinear system. The comparison results 
of closed-loop system performance with synthesized LPV controllers correspond to the analysis results of 
quasi-LPV models. 

Povitskv, Alex: Three-dimensional flow in cavity at yav). ICASE Report No. 2001-31 , (NASA/CR-2001- 
211232), October 24, 2001, 23 pages. AIAA Paper No. 01-2847. 

This study is motivated by three-dimensional flows about protrusions and cavities with an arbitrary 
angle between the external flow and rigid elements. The novel type of a 77 building block 77 cavity flow is 
proposed where the cavity lid moves along its diagonal (Case A). The proposed case is taken as a typical 
representative of essentially three-dimensional highly separated vortical flows having simple single-block 
rectangular geometry of computational domain. Computational results are compared to the previous studies 
where the lid moves parallel to the cavity side walls (Case B). These 3-D lid-driven cavity flows are studied 
by numerical modeling using second-order upwind schemes for convective terms. The volume and plane 
integrals of primary and transversal momentum are introduced to compare cases in a quantitative way. For 
the laminar flow in the cubic cavity, the integral momentum of the secondary flow (which is perpendicular 
to the lid direction) is about an order of magnitude larger than that in Case B. In Case A, the number 
of secondary vortices substantially depends on the Re number. The secondary vortices in the central part 
of the cavity in Case A distinguishes it from Case B, where only corner secondary vortices appear. For a 
rectangular 3-D 3:1:1 cavity the integral momentum of the secondary flow in Case A is an order of magnitude 
larger than that in the benchmark cases. The flowfield in Case A includes a curvilinear separation line and 
non-symmetrical vortices which are discussed in the paper. The estimated Gortler number is approximately 
4.5 times larger in Case A than that in Case B for the same Re number. This indicates that in Case A the 
flow becomes unsteady for smaller Re numbers than in Case B. For developed turbulent flow in the cubic 
cavity, the yaw effect on amplification of secondary flow is as strong as that for the ’aminar flow despite the 
more complex vortical flow pattern in benchmark case B. 
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Parikh, St avail M.: An analysis mechanism for automation in terminal area. ICASE Report No. 2001-32. 
( N AS A /CR- 200 1-2 11 235), October 31, 2001, 15 pages. 

This paper describes an Air Traffic Simulator (ATS) that provides a simulation capability for terminal 
area automation. The ATS is capable of realistic terminal area modeling, automation performance and 
feasibility studies, and system wide fault-tolerance analysis. The ATS is evaluated through modeling of 
self-spacing and self-merging of aircraft following Standard Terminal Arrival (STARs). These algorithms are 
described in detail and preliminary analysis results are provided. The analysis exposes the limitations of 
the algorithms due to their integration. Possible modifications to overcome these limitations are suggested. 
Future work planned for the ATS is described. Overall the ATS is an easily extensible and powerful tool for 
preliminary analysis of new technologies. 

Shin, Jong-\eob. Gary J. Balas, and Alpay M. Kaya: Blending methodology of linear parameter varying 
control synthesis of F-16 aircraft systems. ICASE Report No. 2001-33 , (NASA/CR-2001-211237), October 
29. 2001, 26 pages. To be submitted to the Journal of Guidance, Control, and Dynamics. 

This paper presents the design of a linear parameter varying (LPV) controller for the F-1G longitudinal 
axes over the entire flight envelope using a blending methodology which lets an LPV controller preserve 
performance level over each parameter subspace and reduces computational costs for synthesizing an LPV 
controller. Three blending LPV controller synthesis methodologies are applied to control F-16 longitudinal 
axes. Using a function substitution method, a quasi-LPV model of the F-16 longitudinal axes is constructed 
from the nonlinear equations of motion over the entire flight envelope, including non-trim regions, to facilitate 
synthesis of LPV controllers for the F-16 aircraft. The nonlinear simulations of the blending LPV controller 
show that the desired performance and robustness objectives are achieved across all altitude variations. 

Prieto, Manuel, Ruben S. Montero, and Ignacio M. Llorente: A parallel multigrid solver for viscous flows on 
anisotropic structured grids. ICASE Report No. 2001-34 , (NASA/CR-2001-211238), October 31, 2001. 24 
pages. To be submitted to Parallel Computing. 

This paper presents an efficient parallel multigrid solver for speeding up the computation of a 3-D model 
that treats the flow of a viscous fluid over a flat plate. The main interest of this simulation lies in exhibiting 
some basic difficulties that prevent optimal multigrid efficiencies from being achieved. As the computing 
platform, we have used Coral, a Beowulf-class system based on Intel Pentium processors and equipped with 
GigaNet cLAN and switched Fast Ethernet networks. Our study not only examines the scalability of the 
solver but also includes a performance evaluation of Coral where the investigated solver has been used to 
compare several of its design choices, namely, the interconnection network (GigaNet versus switched Fast- 
Ethernet) and the node configuration (dual nodes versus single nodes). As a reference, the performance 
results have been compared with those obtained with the NAS-MG benchmark. 

Huyse, Luc, and Robert W . Walters: Random field solutions including boundary condition uncertainty for the 
steady-state generalized Burgers equation. ICASE Report No. 2001-35 , (NAS A/CR-2001-21 1239), October 
29, 2001, 30 pages. To be submitted to the Journal of Fluid Mechanics. 

CFD results are subject to considerable uncertainty associated with the operating conditions. Even 
when the operational uncertainty is omitted under very controlled circumstances during wind tunnel exper- 
iments, substantial disagreement between experimental and CFD results persists. This discrepancy must be 
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attributed to model uncertainty. This report discusses the various sources of uncertainty. The need for ad- 
vanced uncertainty modeling is illustrated by means of a computationally inexpensive 1-D Burgers equation 
model. We specifically address the uncertainty due to missing variables (inexact or incomplete differential 
equations). To this extent a random field model is used for the viscosity and the fundamental differences 
between the solutions of the stochastic differential equations and a simple random variable model is high- 
lighted. The Burgers equation theoretically needs to be integrated over an infinite domain. In a deterministic 
approach, the integration domain is cut off at some far field boundary. This truncation effectively ignores 
all variability outside this far field boundary. We present a practical treatment for the uncertainty on the 
boundary conditions. The results indicate that ignoring the boundary condition uncertainty dramatically 
underestimates the variance of the velocity in the interior of the domain. 

Little, Mark E., and Martin E. Kordesch: Band-gap engineering in sputter deposited amorphous /microcry stal- 
line Sc x Ga } - x N. ICASE Report No. 2001-36 , (NASA/CR-2001-211241), November 29. 2001, 11 pages. To 
be submitted to the Proceedings of the MRS Fall Meeting. 

Reactive sputtering was used to grow thin films of Sc x Ga\- x N with scandium concentrations of 20%-70% 
on quartz substrates at temperatures of 300-675 K. X-ray diffraction (XRD) of the films showed either weak 
or no structure, suggesting the films are amorphous or rnicrocrystalline. Optical absorption spectra were 
taken of each sample and the optical band gap was determined. The band gap varied linearly with increasing 
Ga concentration between 2.0 and 3.5 eV. Ellipsometry was used to confirm the band gap measurements and 
provide optical constants in the range 250-1200 nm. ScN and GaN have different crystal structures (rocksalt 
and wurzite, respectively), and thus may form a heterogeneous mixture as opposed to an alloy. Since the 
XRD data were inconclusive, bilayers of ScN/GaN were grown and optical absorption spectra taken. A 
fundamental difference in the spectra between the bilayer films and alloy films was seen, suggesting the films 
are alloys, not heterogeneous mixtures. 

Hart, Roger C., G.C. Herring, and R. Jeffrey Balia: Common-path heterodyne laser-induced thermal acoustics 
for seedless laser velocimctry. ICASE Report No. 2001-37 , (NASA/CR-2001-211252), March 7, 2002, 12 
pages. To be submitted to Optic Letters. 

We demonstrate the use of a novel technique for the detection of heterodyne laser-induced thermal 
acoustics signals, which allows the construction of a highly stable seedless laser velocimeter. A common-path 
configuration is combined with quadrature detection to provide flow' direction, greatly improve robustness 
to misalignment and vibration, and give reliable velocity measurement at low r flow* velocities. Comparison 
with Pitot tube measurements in the freestream of a windtunnel shows root-mean-square errors of 0.67 m/s 
over the velocity range 0-55 m/s. 

Zhang, Yong-Tao, and Chi- Wang Shu: High order WENO schemes for Hamilton- Jacobi equations on trian- 
gular meshes . ICASE Report No. 2001-38 , (NASA/CR-2001-211256), December 26, 2001, 36 pages. Sub- 
mitted to the SIAM Journal on Scientific Computing. 

In this paper we construct high order weighted essentially non-oscillatorv (WENO) schemes for solving 
the nonlinear Hamilton- Jacobi equations on two-dimensional unstructured meshes. The main ideas are nodal 
based approximations, the usage of monotone Hamiltonians as building blocks on unstructured meshes, non- 
linear weights using smooth indicators of second and higher derivatives, and a strategy to choose diversified 
smaller stencils to make up the bigger stencil in the WENO procedure. Both third-order and fourth-order 
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WENO schemes using combinations of second-order approximations with nonlinear weights are constructed. 
Extensive numerical experiments are performed to demonstrate the stability and accuracy of the meth- 
ods. High-order accuracy in smooth regions, good resolution of derivative singularities, and convergence to 
viscosity solutions are observed. 

Xu, Kun: Regularization of the Chapman- Enskog expansion and its description of shock structure . ICASE 
Report No. 2001-39 , (NASA/CR-2001-211268), December 26, 2001, 11 pages. To be submitted to Physics 
of Fluids. 

In the continuum transition flow regime, we propose to truncate the Chapman-Enskog expansion of the 
Boltzmann equation to the Navier-Stokes order only without going to the Burnett or super Burnett orders. 
However, the local particle collision time has to be generalized to depend not only on the local macroscopic 
flow variables, but also their gradients in the rarefied gas regime. Based on the gas-kinetic BGK model, 
the relation between the conventional collision time and the general one is obtained. More specificially. 
a generalized constitutive relation for stress and heat flux is proposed. This new model is applied to the* 
study of argon gas shock structure. There is good agreement between the predicted shock structure and 
experimental results for a wide range of Mach numbers. 

Hu, Fang Q., and Harold L. Atkins: Eigensolution analysis of the discontinuous Galerkin method with 
non-uniform grids , part I: One space dimension . ICASE Report No. 2001-40 , (NASA/CR-2001-211269). 
December 26, 2001, 35 pages. To be submitted to the Journal of Computational Physics. 

We present a detailed study of spatially propagating waves in a discontinuous Galerkin scheme applied to 
a system of linear hyperbolic equations. We start with an eigensolution analysis of the semi-discrete system 
in one space dimension with uniform grids. It is found that, for any given order of the basis functions, 
there are at most two spatially propagating numerical wave modes for each physical wave of the Partial 
Differential Equations (PDE). One of the modes can accurately represent the physical wave of the PDE and 
the other is spurious. The directions of propagation of these two numerical modes are opposite, and, in most 
practical cases, the spurious mode has a large damping rate. Furthermore, when an exact characteristics 
split flux formula is used, the spurious mode becomes non-existent. For the physically accurate mode, it is 
shown analytically that the numerical dispersion relation is accurate to order 2p 4* 2 where p is the highest 
order of the basis polynomials. The results of eigensolution analysis are then utilized to study the effects 
of a grid discontinuity, caused by an abrupt change in grid size, on the numerical solutions at either side of 
the interface. It is shown that, due to "mode decoupling," numerical reflections at grid discontinuity, when 
they occur, are always in the form of the spurious non-physical mode. Closed form numerical reflection and 
transmission coefficients are given and analyzed. Numerical examples that illustrate the analytical findings 
of the paper are also presented. 

Jow, Kenny G., and Theo J. Dingemans: Liquid crystals derived from 2-phenyl-isoindoles. ICASE Report No. 
2QQ1-4E (NASA/CR-2001-211270), January 15, 2002, 18 pages. To appear in Liquid Crystals. 

2-Phenyl-isoindole was investigated as the rigid core unit in a series of asymmetric mesogenic molecules. 
When the 2-phenyl-isoindole core was terminated with a hexyl tail, no mesophase formation could be ob- 
served. When 4-n-(tridecafluorohexyl) was used, however, we observed both monotropic and enantiotropic 
phase behavior. We found that most functionalities at the anhydride 5-position results in the formation of 
smectic* A (SmA) phases in the temperature range of 70-180°C. Functionalities at the anhydride 4-position 
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suppress mesophase formation. Large substituents (-Br, -NO-?) and symmetric substitution patterns (5,6- 
dichloro, 4,7-didiloro and 4,5,6,7-tetrachloro) on the anhydride moiety increase the melting point and desta- 
bilize the mesophase. Temperature dependent X-ray diffraction experiments suggest an interdigitated SmA 
packing for this family of compounds. 

Heggernes. P., S.C. Eisenstat, G. Kumfcrt, and A. Pothen: The computational complexity of the Mini- 
mum Degree algorithm. ICASE Report No. 2001-42 , (NASA/CR-2001-211421), January 15. 2002. 16 pages. 
Proceedings of the Norwegian Computer Science Conference NIK 2001. 

The Minimum Degree algorithm, one of the classical algorithms of sparse matrix computations, is widely 
used to order graphs to reduce the work and storage needed to solve sparse systems of linear equations. There 
has been extensive research involving practical implementations of this algorithm over the past two decades. 
However, little has been done to establish theoretical bounds on the computational complexity of these 
implementations. We study the Minimum Degree algorithm, and prove time complexity bounds for its 
widely used variants. 

Harrison, J.S., and Z. Ounaies: Piezoelectric polymers . ICASE Report No. 2001-43 , (NASA/CR-2001- 
211422), March 28, 2002, 31 pages. To appear in the Encyclopedia of Smart Materials, John Wiley, December 
2001. 

The purpose of this review is to detail the current theoretical understanding of the origin of piezoelec- 
tric and ferroelectric phenomena in polymers: to present the state-of-the-art in piezoelectric polymers and 
emerging material systems that exhibit promising properties; and to discuss key characterization methods, 
fundamental modeling approaches, and applications of piezoelectric polymers. Piezoelectric polymers have 
been known to exist for more than forty years, but in recent years they have gained notoriety as a valuable 
class of smart materials. 

Carpenter, Mark H., David Gottlieb, and Chi- Wang Shu: On the conservation and convergence to weak 
solutions of global schemes . ICASE Report No. 2001-44 , (NASA/CR-2001-211423), March 1. 2002, 20 pages. 
To appear in the Journal of Scientific Computing. 

In this paper we discuss the issue of conservation and convergence to weak solutions of several global 
schemes, including the commonly used compact schemes and spectral collocation schemes, for solving hy- 
perbolic conservation laws. It is shown that such schemes, if convergent boundedly almost everywhere, will 
converge to weak solutions. The results are extensions of the classical Lax-Wendroff theorem concerning 
conservative schemes. 

Wu, Xuesong, and Meelan Choudhari: Linear and nonlinear instabilities of Blasius boundary layer per- 
turbed by streamwise vortices part II: Intermittent instability induced by long-wavelength Klebanoff modes. 
ICASE Report No. 2001-45 , (NASA/CR-2001-211424), March 8, 2002, 46 pages. To be submitted to the 
Journal of Fluid Mechanics. 

We present theoretical results on the stability properties of a Blasius boundary layer perturbed by 
Klebanoff distortions with a relatively long spanwise scale. Even relatively weak Klebanoff modes can alter 
the near- wall curvature of the underlying flow by 0(1) and, hence, introduce linear instabilities with larger 
characteristic growth rates and frequencies than those of the Tollmien-Schlichting waves in an unpertubed 
Blasius flow. A localised distortion supports both sinuous and varicose modes of instability, with the growth 
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rates of the sinuous modes being likely to be larger, in general. Overall, the instability is intermittent in 
time and localised in space, being confined to a finite part of the Klebanoff mode cycle and to a specific 
window (s) along the streamwise direction. A spanwise periodic distortion supports spatially quasi-periodie 
modes (via the parametric resonance mechanism), which may be viewed as modified T-S waves with excess 
growth rates when the Klebanoff modes are weak. In spite of the simplifications involved in this theory, its 
predictions appear qualitatively consistent with some of the unusual characteristics of the high-frequency 
wavepackets observed during previous experiments. The nonlinear development of a localised sinuous mode 
is followed across a sequence of asymptotic regimes using the non-equlibrium critical-layer theory. 

Walters, Robert W., and Luc Huyse: Uncertainty analysis for fluid mechanics with applications. ICASE 
Report No. 2002-1 , (NASA/CR-2002-211449), March 29, 2002, 50 pages. 

This paper reviews uncertainty analysis methods and their application to fundamental problems in fluid 
dynamics. Probabilitic (Monte-Carlo, Moment methods, Polynomial Chaos) and non-prohabilistic meth- 
ods (Interval Analysis, Propagation of error using sensitivity derivatives) are described and implemented. 
Results are presented for a model convection equation with a source term, a model non-linear convection- 
diffusion equation; supersonic flow over wedges, expansion corners, and an airfoil; and two-dimensional 
laminar boundary layer flow. 

Rubinstein, R. , and Guowei He: The Eulerian time correlation function in homogeneous isotropic turbulence. 
ICASE Report No. 2002-2 , (NASA/CR-2002-211450), March 29, 2002, 14 pages. To appear in the Proceed- 
ings of the International Workshop on Statistical Theories and Computational Approaches to Turbulence. 

Two general models are proposed for the Eulerian time correlation function in homogeneous isotropic 
turbulence. The first is based on continued fraction approximations to its Laplace transform, and the second 
is based on random sweeping by a possibly non-Gaussian velocity field. Both models can give reasonable 
quantitative agreement with DNS data for moderate time separations over which the time correlation func- 
tions at different wavenumbers exhibit a common self-similar form. 

Geser, Alfons: Loops of superexponential lengths in one-rule string rewriting. ICASE Report No. 2002-3. 
(NASA/ CR-2002-21 1451) , March 27, 2002, 16 pages. To be submitted to the International Conference on 
Rewriting Techniques and Applications. 

Loops are the most frequent cause of non-termination in string rewriting. In the general case, non- 
terminating, non-looping string rewriting systems exist, and the uniform termination problem is undecidable. 
For rewriting with only one string rewriting rule, it is unknown whether non-terminating, non-looping systems 
exist and whether uniform termination is decidable. If in the one-rule case, non-termination is equivalent 
to the existence of loops, as McNaughton conjectures, then a decision procedure for the existence of loops 
also solves the uniform termination problem. As the existence of loops of bounded lengths is decidable, the 
question is raised how long shortest loops may be. We show that string rewriting rules exist whose shortest 
loops have superexponential lengths in the size of the rule. 
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Krueger, Ronald, and Pierre J. Minguet: Influence of 2D finite element modeling assumptions on debonding 
prediction for composite skin-stiffener specimens subjected to tension and bending. ICASE Report No. 2002-4 , 
(NASA /CR..-2002-21 1452) , March 29, 2002, 35 pages. To be presented at the International Conference on 
Composite Structures. 

The influence of two-dimensional finite element modeling assumptions on the debonding prediction 
for skin-stiffener specimens was investigated. Geometrically nonlinear finite element analyses using two- 
dimensional plane-stress and plane-strain elements as well as three different generalized plane strain type 
approaches were performed. The computed deflections, skin and flange strains, transverse tensile stresses 
and energy release rates were compared to results obtained from three-dimensional simulations. The study 
showed that for strains and energy release rate computations the generalized plane strain assumptions yielded 
results closest to the full three-dimensional analysis. For computed transverse tensile stresses the plane stress 
assumption gave the best agreement. Based on this study it is recommended that results from plane stress 
and plane strain models be used as upper and lower bounds. The results from generalized plane strain 
models fall between the results obtained from plane stress and plane strain models. Two-dimensional models 
may also be used to qualitatively evaluate the stress distribution in a ply and the variation of energy release 
rates and mixed mode* ratios with delamination length. For more accurate predictions, however, a three- 
dimensional analysis is required. 
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INTERIM REPORTS 


Eidson, Thomas M., and Merrell L. Patrick: ICASE Workshop on Programming Computational Grids. 
ICASE Interim Report No. 38 , (NASA/CR-2001-211224), October 15. 2001, 21 pages. 

A workshop on Programming Computational Grids for distributed applications was held on April 12-13. 
2001 at ICASE, NASA Langley Research Center. The stated objective of the workshop was to define, discuss, 
and clarify issues critical to the advancement of Problem Solving Environments/Computational Frameworks 
for solving large multi-scale, multi-component, scientific applications using distributed, heterogeneous com- 
puting systems. This report documents a set of recommendations for NASA that suggest an approach for 
developing an application development environment that will meet future application needs. 


Munoz, Cesar, and Micaela Mayero: Real automation in the field . ICASE Interim Report No. 39 , (NASA/CR- 
2001-211271), January 22, 2002, 14 pages. 

We provide a package of strategies for automation of non-linear arithmetic in PVS. In particular, we 
describe a simplification procedure for the field of real numbers and a strategy for cancellation of common 
terms. 
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ICASE COLLOQUIA 


Name/ Affiliation/Title 


Achi Brandt. The Weizmann Institute of Science. Israel 
“Multiscale Computational Methods for Molecular Dynamics” 

Ken Schwaher, Advanced Development Methods, Inc. 

ICASE Series on Modern Programming Practices: “Agile Processes: A Case 
Study of a Scrum Project” 

Gregory Carman, University of California, Los Angeles 
ICASE Series on Morphing: “Investigating the Passive Damping Properties 
of Active Materials” 

Michael Hill, Object Mentor, Inc. 

ICASE Series on Modern Programming Practices: “The Four Keys to 
Test-first Design: Laziness, Resentment, Xenophobia, and Guile” 

George Haller, Brown University 
“Lagrangian Coherent Structures and their Control” 

Thomas Bewlev, University of California, San Diego 
ICASE Series on Morphing: “Adjoint and Riccati: Essential Tools in the 
Analysis and Control of Transitional and Turbulent Flow Systems ' 

Terrence Weisshaar, Purdue University 
ICASE Series on Morphing: “Aeroelastic Tailoring for Energy Efficient 
Morphing Aircraft - Finding the Right Stuff’ 

Dieter Wolf, Argonne National Laboratory 
ICASE Series on Nanotechnology: “Plastic Deformation of Nanocrystalline 
Materials by Molecular-dynamics Simulation” 

Inderjit Chopra, University of Maryland, College Park 
ICASE Series on Morphing: “Review on the Status of Application of Smart 
Structures Technology to Rotorcraft Systems” 

Michael Two, Thought Works, Inc. 

ICASE Series on Modern Programming Practices: “Continuous Integration” 


Date 


October 2 


October 2 


October 3 


October 4 


October 22 


October 23 


October 25 


October 30 


October 31 


November 5 
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Name/ Affiliation/Title 


Date 


Joel Grasmeyer, AeroVironment, Inc. 

ICASE Series on Morphing: “Perspectives on Micro Air Vehicles and 
Electronic Packaging” 

Tony Ingraffea, Cornell University 

ICASE Series on Nanotechnology: “A Multiscale Modeling Approach to 
Crack Initiation in Aluminum Polycrystals” 

David Geohegan, Oak Ridge National Laboratory 
“Time-resolved Growth Studies and Applications of Single-wall Carbon 
Nanotubes” 

Richard James, University of Minnesota 
ICASE Series on Nanotechnology: “Deforming Films of Active Materials: 

New Concepts for Producing Motion at Small Scales using Applied Fields” 

Reeti Katoch- Rouse, Lehigh University 
“Controlled Reductive Cydization for the Synthesis of dihydrobenzofurans/ 
Stereocontrolled Preparation of Functionalized Dienes for Decalin Synthesis” 

Deborah Mayo, Virginia Polytechnic Institute and State University 
ICASE Series on Risk-based Design: “Learning from Error and an Error 
Statistical Philosophy of Science” 

Anatole Beck, University of Wisconsin, Madison 
“Some Interesting Anomalies in the Uniqueness of Solutions of Differential 
Equations” 

Armen Der Kiureghian, University of California, Berkeley 
ICASE Series on Risk-based Design: “Algorithms for Reliability-based Optimal 
Structural Design” 

David McDowell, Georgia Institute of Technology 
ICASE Series on Nanotechnology: “Computational Modeling of Plasticity and 
Fracture at Microstructure Scales” 

Mohammad Khalessi, Unipass Technologies 
ICASE Series on Risk-based Design: “Probabilistic Technology at 
PredictionProbe, Inc.” 
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November 9 


November 13 


November 14 


November 15 


November 16 


November 16 


November 20 


November 26 


November 29 


December 3 



Name/ Affiliation/Title 


Date 


Satish Kumar, Georgia Institute of Technology 
“Films, Fibers, and Coatings from Carbon Nano Tubes, Nano Fibers, arid their 
Composites with Polymers" 

Kyung Choi, The University of Iowa 

ICASE Series on Risk-based Design: “Advances in Reliability-based Design 
Optimization and Probability Analysis - Part I” 

Kyung Choi, The University of Iowa 

ICASE Series on Risk-based Design: “Advances in Reliability-based Design 
Optimization and Probability Analysis - Part II” 

Harry Dorn, Virginia Polytechnic Institute and State University 
ICASE Series on Multifunctional Materials: “Nanoscale Carbonaceous Building 
Blocks: Properties and Applications” 

Do Young, Kwak KAIST, Korea 

“Convergence of Multigrid for Nonconforming FE and Cell-centered FD” 

Kie Eom, The George Washington University 
“Random Field Models for Analysis of Natural and Artificial Patterns” 

Emmanuel Giannelis, Cornell University 
“New Advances in Polymer Nanocomposites” 

Boris Berkovski, Global Science and Technology, Inc. 

“Some New Approaches in Aeronautical Modeling” 

Michael Trosset, The College of William and Mary 
“Uncertainty in Computational Science: First Impressions” 

Brien Alkire, University of California, Los Angeles 
“Convex Optimization Problems Involving Autocorrelation Sequences” 

Roger Ghanem, The Johns Hopkins University 
ICASE Series on Risk-based Design: “Propagation and Management of 
Uncertainty in Mechanics-based Models” 

Nicholas Bailey, Cornell University 

ICASE Series on Nanotechnology: “Aspects of Mixed Atomistic-continuum 
Modeling” 


December 5 


December 10 


December 11 


December 13 


January 18 


January 25 


February 5 


February 6 


February 8 


February 11 


February 14 


February 22 
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Name/Affiliation/Title 


Date 


Dmitry Keselman, Atraxis North America. Inc. February 25 

“On Some Problems of Scheduling and the Airline Network Analysis' 5 

Luis Crespo, University of Delaware March 4 

“Deterministic and Stochastic Optimal Control via Bellman’s Principle 51 

John Gilbert, Palo Alto Research Center and MIT Laboratory for Computer Science March 7 
“Smart Matter: Frontiers of Computation 55 

George Lauder, Harvard University March 15 

ICASE Series on Morphing: “Fish Morphing: Experimental Hydrodynamics of 
Locomotion in Fishes 55 

Larry Howell, Brigham Young University March 25 

“Compliant Micro Mechanism Analysis and Design 55 

John Renaud, University of Notre Dame March 28 

ICASE Series on Risk-based Design: “Decision-based Collaborative Optimization 
and Estimation of Propagated Uncertainty 55 
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OTHER ACTIVITIES 


On November 6 7. 2001, ICASE/USRA and NASA Langley co-sponsored a workshop on Revolutionary 
Aerospace Systems Concepts for Human/Robotic Exploration of the Solar System at the Quality Inn and 
the Holiday Inn in Hampton, Virginia. The primary focus of this workshop was to understand how, with the 
incorporation of revolutionary aerospace systems concepts over the next 10 40 years, humans and machines 
can be synergistically combined to physically and virtually reduce the time and distance barriers associated 
with exploring beyond low Earth orbit. There were 85 attendees and a NASA Conference Proceedings will 
be published. 
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ICASE STAFF 


I. ADMINISTRATIVE 

Manuel D. Salas, Director, M.S., Aeronautics and Astronautics, Polytechnic Institute of Brooklyn, 1970. 
Fluid Mechanics and Numerical Analysis. 

Linda T. Johnson, Office and Financial Administrator 

Barbara A. Cardasis, Administrative Secretary 

Etta M. Morgan, Accounting Supervisor 

Emily N. Todd, Conference Manager /Executive Assistant 

Shannon K. Verstynen, Information Technologist 

Gwendolyn W. Wesson, Contract Accounting Clerk 

Shouben Zhou, Systems Manager 

II. SCIENCE COUNCIL 

David Gottlieb, (Chair) Professor. Division of Applied Mathematics, Brown University. 

II han Aksay, Professor, Engineering Quad, Princeton University. 

Lee Beach, Professor, Department of Physics, Computer Science & Engineering, Christopher Newport Uni- 
versity. 

Jack Dongarra, Distinguished Professor. Department of Computer Science, University of Tennessee. 

Joseph E. Flaherty, Amos Eaton Professor, Departments of Computer Science and Mathematical Sciences. 
Rensselaer Polytechnic Institute. 

Forrester Johnson, Aerodynamics Research, Boeing Commercial Airplane Group. 

John C. Knight, Professor, Department of Computer Science, School of Engineering and Applied Science. 
University of Virginia. 

Robert W. MacCormack, Professor, Department of Aeronautics and Astronautics, Stanford University. 

Stanley G. Rubin, Professor. Department of Aerospace Engineering and Engineering Mechanics, University 
of Cincinnati. 

Manuel D. Salas, Director, ICASE, NASA Langley Research Center. 

46 


i 



III. RESEARCH FELLOWS 


Dimitri Mavriplis - Ph.D., Mechanical and Aerospace Engineering, Princeton University, 1988. Applied k 
Numerical Mathematics [Grid Techniques for Computational Fluid Dynamics]. (February 1997 to August 
2005) 

.losip Loncaric- Ph.D., Applied Mathematics, Harvard University, 1985. Applied k Numerical Mathematics 
[Multidisciplinary Design Optimization]. (March 2001 to August 2002) 


IV. SENIOR STAFF SCIENTISTS 


Brian G. Allan - Ph.D., Mechanical Engineering, University of California at Berkeley, 1996. Applied k 
Numerical Mathematics [Multidisciplinary Design Optimization], (February 1996 to November 2003) 

Maria Consiglio - M.S., Computer Science, University of Illinois-Urbana, 1982. Computer Science [Crew 
Systems Research for Aviation Capacity and Safety]. (June 2001 to June 2003) 

Thomas M. Eidson - Ph.D., Mechanical Engineering, University of Michigan, 1982. Computer Science 
[Distributed Computing]. (October 2000 to September 2003) 

Alfons E. Geser - Ph.D., Computer Science, University of Passau, Germany, 1991. Computer Science [Formal 
Methods], (January 2001 to December 2002) 

Vasyl M. Harik - Ph.D., Physics, Moscow State University, 1993 and Mechanical Engineering, University 
of Delaware, 1997. Structures k Materials [Composites and Failure Mechanics]. (October 2000 to October 

2003) 

Roger C. Hart - Ph.D., Physics, University of Tennessee, 1991. Fluid Mechanics [Measurement Science and 
Technology]. (December 1998 to October 2003) 

Guowei He - Ph.D., Theoretical and Applied Mechanics, Northwestern Polytechnic University, Xian, China, 
1991. Fluid Mechanics [Turbulence Modeling and Direct Numerical Simulation]. (July 2000 to June 2003) 

Navin Jaunky - Ph.D., Mechanical Engineering, Old Dominion University, 1995. Structures k Materials 
[Composite Structural Damage Tolerance and Residual Strength Methodologies]. (January 2001 to December 

2004) 

Li-Shi Luo - Ph.D., Physics, Georgia Institute of Technology, 1993. Computer Science [Parallel Algorithms], 
(November 1996 to August 2004) 

Zoubeida Ounaies - Ph.D., Engineering Science and Mechanics, The Pennsylvania State University, 1995. 
Structures k Materials [Characterization of Advanced Piezoelectric Materials], (March 1999 to December 
2001 ) 

Cheol Park - Ph.D., Macromolecular Science and Engineering, The University of Michigan, 1997. Structures 
k Materials [Electro-active Materials], (November 2000 to October 2002) 
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Alexander Povitskv - Ph.D., Mechanical Engineering, Moscow Institute of Steel and Alloys Technology 
(MISA), Russia, 1988. Computer Science [Parallelization and Formulation of Higher Order Schemes for 
Aeroacoustics Noise Propagation]. (October 1997 to November 2001) 


V. SCIENTIFIC STAFF 

Scott C. Beeler - Ph.D., Applied Mathematics, North Carolina State University, 2000. Applied & Numerical 
Mathematics [Nonlinear Suboptimal Feedback Control]. (March 2001 to February 2003) 

Fred E. Blough - B.A., Business Administration, Upper Iowa University, 1980. Business Development 
Specialist [Stimulation and Acceleration of New Technology and Intellectual Capital for LaRC]. (February 
2002 to February 2004) 

Frank Bussink - B.Sc., Aeronautical Engineering/Computer Science, Polytechnics of Amsterdam, The Nether- 
lands, 1996. Computer Science [Distributed Air-Ground Traffic Management,]. (January 2002 to December 
2003) 

Theodorus Dingemans - Ph.D., Organic Chemistry, University of North Carolina at Chapel Hill, 1998. 
Structures & Materials. (September 2000 to August 2002) 

Boris Diskin - Ph.D., Applied Mathematics, The Weizmann Institute of Science, Israel, 1998. Applied k 
Numerical Mathematics [Convergence Acceleration]. (July 1998 to September 2004) 

Alicia M. Dwyer - M.S., Mechanical Engineering, Aerospace Track, The George Washington University, 2001. 
Applied & Numerical Mathematics [Planetary Exploration]. (July 2001 to July 2003) 

Sarah- Jane V. Frankland - Ph.D., Chemistry (Physical), The Pennsylvania State University, 1997. Structures 
& Materials [Computational Nanotechnology], (June 2001 to June 2003) 

Hanne Gottliebsen - Ph.D., Computer Science, University of St. Andrews, Scotland, 2001. Computer Science 
[Formal Methods]. (August 2001 to August 2003) 

Jill L. Hanna - M.S., Aerospace Engineering (Astronautics), The George Washington University, 2001. 
Applied & Numerical Mathematics [Planetary Exploration]. (July 2001 to July 2003) 

Luc Huyse - Ph.D., Civil Engineering, Structures, University of Calgary, Canada, 1999. Applied &: Numerical 
Mathematics [Managing Uncertainties in Multidisciplinary Research]. (October 1999 to November 2001) 

Ronald Krueger - Ph.D., Aerospace Engineering, University of Stuttgart, Germany, 1996. Structures & 
Materials [Analysis of Composite Delamination of Structures]. (August 2000 to August 2003) 

Mark E. Little - Ph.D., Physics, Ohio University, 2001. Structures & Materials [Intelligent Optics]. (April 
2001 to April 2003) 

Cesar A. Munoz - Ph.D., Computer Science, University of Paris 7, 1997. Computer Science [Formal Methods 
Research for Safety Critical Systems]. (May 1999 to April 2003) 
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Gregory M. Odegard - Ph.D., Engineering, University of Denver, 2000. Structures k Materials [Computa- 
tional Nanotechnology]. (February 2002 to January 2004) 

Jason H. Rouse - B.S., Chemistry, Lehigh University, 1995. Structures k Materials [Molecular Self-assembly]. 
(May 2001 to May 2003) 

Jong-Yeob Shin - Ph.D., Aerospace Engineering, University of Minnesota, 2000. Research Assistant, Uni- 
versity of Minnesota. Applied k Numerical Mathematics [Advanced Control Methods]. (November 2000 to 
October 2002) 

Kent A. Watson - Ph.D., Organic/ Polymer Chemistry, Virginia Commonwealth University, 1998. Structures 
k Materials [Nanocomposites]. (April 2001 to March 2003) 

David W. Way - Ph.D., Aerospace Engineering, Georgia Institute of Technology, 2001. Applied k Numerical 
Mathematics [Planetary Exploration]. (August 2001 to August 2003) 

Kristopher Wise - Ph.D., Theoretical/ Computational Chemistry, University of Oklahoma, 1999. Structures 
k Materials [Computational Chemistry for Nanostructured Materials]. (November 2001 to October 2003) 

Tian-Bing Xu - Ph.D., Materials Engineering, The Pennsylvania State University, 2002. Structures k 
Materials [Electro Active Materials], (January 2002 to January 2004) 

VI. VISITING SCIENTISTS 

Kab Seok Kang - Ph.D., Mathematics, Korea Advanced Institute of Science k Technology, 1999. Post- 
doctoral Research Scientist, Korea Advanced Institute of Science k Technology. Applied k Numerical 
Mathematics [Multigrid Algorithms for Partial Differential Equations Discretized on Unstructured Grids]. 
(February 2001 to December 2002) 


VII. SHORT-TERM VISITING SCIENTISTS 

Saul Abarbanel - Ph.D., Theoretical Aerodynamics, Massachusetts Institute of Technology, 1959. Professor, 
Department of Applied Mathematics, Tel Aviv University, Israel. Applied k Numerical Mathematics [Global 
Boundary Conditions for Aerodynamics and Aeroacoustic Computations]. (October 2001) 

Jaehwan Kim - Ph.D., Engineering Science and Mechanics, The Pennsylvania State U niversity, 1995. As- 
sistant Professor, Department of Mechanical Engineering, Inha University, Incheon, Korea. Structures k 
Materials [Smart Materials and Flow Control]. (December 2001 to February 2002) 

Hans Peter Monner - Ph.D.. Adaptronics, Technical University of Braunschweig, Germany, 2000. Head, 
Department of Adaptronics, Institute of Structural Mechanics of the German Aerospace Center. Structures 
k Materials. (March 2002 to June 2002) 

Sun Mok Paik - Ph.D., Physics, University of Maryland, 1988. Assistant Associate Professor, Department of 
Physics, Kangwon National University, Korea. Structures k Materials [Computational Materials]. (January 
2002 to February 2002) 
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Avi Seifert - Ph.D., Fluid Mechanics, Tel- Aviv University, 1990. Senior Lecturer, Department of Fluid 
Mechanics k Heat Transfer, Faculty of Engineering, Tel Aviv University. Fluid Mechanics. (February 2002) 


Robert Walters - Ph.D.. Aerospace Engineering, North Carolina State University, 1984. Research Profes- 
sor, Department of Aerospace and Ocean Engineering, Virginia Polytechnic Institute and State University. 
Applied k Numerical Mathematics [Managing Uncertainties in Multidisciplinary Research]. (February 2002) 


VIII. ASSOCIATE RESEARCH FELLOW 

David E. Keyes - Ph.D., Applied Mathematics, Harvard University, 1984. Computer Science [Parallel Nu- 
merical Algorithms] 


IX. CONSULTANTS 

Tomasz Arciszewski - Ph.D., Technical Sciences, Warsaw University of Technology, Poland. 1975. Profes- 
sor, Civil, Environmental, and Infrastucture Program, School of Information Technology and Engineering, 
George Mason University. Applied k Numerical Mathematics [Evolutionary Computation and Complexity 
in Innovative Design] 

H. Thomas Banks - Ph.D., Applied Mathematics, Purdue University, 1967. Professor, Department of Math- 
ematics, Center for Research in Scientific Computations, North Carolina State University. Applied k Nu- 
merical Mathematics [Control Theory] 

Przemyslaw Bogacki - Ph.D., Mathematical Sciences, Southern Methodist University, 1990. Assistant Profes- 
sor, Department of Mathematics and Statistics, Old Dominion University. Applied k Numerical Mathematics 
[High Performance Methods in Nontraditional CFD] 

Achi Brandt - Ph.D., Mathematics, The Weizmann Institute of Science, 1965. Professor, Department of 
Applied Mathematics, The Weizmann Institute of Science, Israel. Applied k Numerical Mathematics [Con- 
vergence Acceleration] 

Timothy D. Bryant - No College Degree. TDB Engineering, Gloucester, VA. Structures k Materials [Design 
and Fabrication of Standard Configuration for Making and Testing Sensor and/or Actuator Elements from 
Electroactive Polymers] 

Thomas \\ . Crockett - B.S., Mathematics, The College of William k Mary, 1977. Senior Research Associate, 
Computational Science Cluster, The College of William k Mary. Computer Science [Scientific Visualization] 

Ayodeji O. Demuren - Ph.D., Mechanical Engineering, Imperial College London, United Kingdom, 1979. 
Associate Professor, Department of Mechanical Engineering and Mechanics, Old Dominion University. Fluid 
Mechanics [Numerical Modeling of Turbulent Flows] 

Dave E. Eckhardt - Ph.D., Computer Science, George Washington University, 1978. Retired. Computer 
Science [Operational Concepts of National Aerospace System Needs] 
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Isaac Elishakoff- Ph.D., Mechanical Engineering, Moscow Power Engineering Institute and State University, 
1971. Professor, Department of Mechanical Engineering, Florida Atlantic University. Structures k Materials 
[Reliability-based Structural Design Technology] 

Sharath Girimaji - Ph.D., Mechanical Engineering, Cornell University, 1990. Associate Professor, Depart- 
ment of Aerospace Engineering, Texas A&M University. Fluid Mechanics [Turbulence and Combustion] 

David Gottlieb - Ph.D., Numerical Analysis, Tel-Aviv University, Israel, 1972. Ford Foundation Professor k 
Chair, Division of Applied Mathematics, Brown University. Applied k Numerical Mathematics [Boundary 
Conditions for Hyperbolic Systems] 

Jan S. Hesthaven - Ph.D., Applied Mathematics/Numerical Analysis, Technical University of Denmark, 1995. 
Visiting Assistant Professor, Division of Applied Mathematics, Brown University. Applied k Numerical 
Mathematics [Computational Electromagnetics] 

Fang Q. Hu - Ph.D., Applied Mathematics, Florida State University, 1990. Assistant Professor, Department 
of Mathematics and Statistics, Old Dominion University. Fluid Mechanics [Aeroacoustics] 

Mohammad R. Khalessi - Ph.D., Mechanics and Structures, University of California-Los Angeles, 1983. 
Chief Product Development Officer, PredictionProbe, Newport Beach, CA. Practical Implementation of 
Probabilistic Technology 

Frank Kozusko- Ph.D., Computational and Applied Mathematics, Old Dominion University, 1995. Assistant 
Professor, Department of Mathematics, Hampton University. Fluid Mechanics [Airfoil Design] 

R. Michael Lewis - Ph.D., Mathematical Sciences, Rice University, 1989. Assistant Professor, Department of 
Applied Mathematics, The College of William k Mary. Applied k Numerical Mathematics [Multidisciplinary 
Optimization and Managing Uncertainties] 

Wu Li - Ph.D., Mathematics, The Pennsylvania State University, 1990. Associate Professor, Department of 
Mathematics and Statistics, Old Dominion University. Applied k Numerical Mathematics [Optimization] 

Hong Zong Lin - Ph.D., Structures, University of California-Berkeley, 1990. Chief Technology Officer, 
PredictionProbe, Newport Beach, CA. Practical Implementaion of Probabilistic Technology 

Gerald Liittgen - Ph.D., Computer Science, University of Passau, Germany, 1998. Professor. Department of 
Computer Science, The University of Sheffield, United Kingdom. Computer Science [Formal Methods] 

Frank T. Lynch - B.S., Aero Engineering, University of Notre Dame, 1955. Lynch Aerodyn Consulting, 
Yorba Linda, CA. Fluid Mechanics [Reynolds Number Scaling Experiences and Lessons Learned for Subsonic 
Transport Aircraft] 

James E. Martin - Ph.D., Applied Mathematics, Brown University, 1991. Assistant Professor, Department 
of Mathematics, Christopher Newport University. Fluid Mechanics [Turbulence and Computation] 

William H. Mason - Ph.D., Aerospace Engineering, Virginia Polytechnic Institute, 1975. Professor, Depart- 
ment of Aerospace and Ocean Engineering and Multidisciplinary Analysis and Design Center for Advanced 
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Vehicles, Virginia Polytechnic Institute and State University. Fluid Mechanics [Development of Roadmap 
Future Research in Aerodynamics] 


Karla Mossi - Ph.D., Mechanical Engineering, Old Dominion University, 1998. Assistant Professor. Depart- 
ment of Mechanical Engineering, Virginia Commonwealth University. Structures & Materials [Electro Active 
Materials] 

Robert S. Nerem - Ph.D., Aerospace Engineering, The University of Texas-Austin, 1989. Associate Professor. 
Department of Aerospace Engineering Sciences, University of Colorado. Revolutionary Systems Concepts 
[Comet /Asteroid Protection System] 

Zoubeida Ounaies - Ph.D., Engineering Science and Mechanics, The Pennsylvania State University, 1995. 
Assistant Professor, Department of Mechanical Engineering, Virginia Commonwealth University. Structures 
& Materials [Electro Active Materials] 

Merrell L. Patrick - Ph.D., Mathematics, Carnegie Institute of Technology. 1964. Retired. Computer Science 

Alex Pothen - Ph.D., Applied Mathematics, Cornell University, 1984. Professor, Department of Computer 
Science, Old Dominion University. Computer Science [Parallel Numerical Algorithms] 

Chi- Wang Shu - Ph.D.. Mathematics, University of California-Los Angeles, 1986. Professor. Division of 
Applied Mathematics, Brown University. Fluid Mechanics [Computational Aeroacoustics] 

Ralph C. Smith - Ph.D., Mathematics, Montana State University, 1990. Professor. Department of Mathe- 
matics, North Carolina State University. Applied &: Numerical Mathematics [Optimal Control Techniques 
for Structural Acoustics Problems] 

Eduardo A. Socolovsky - Ph.D., Mathematics, Carnegie-Mellon University, 1984. Associate Professor, De- 
partment of Mathematics, Hampton University. Applied & Numerical Mathematics [Numerical Methods for 
Visualization and Multimedia] 

Michael W. Trosset. - Ph.D., Statistics, University of California-Berkeley, 1983. Assistant Professor, De- 
partment of Mathematics, The College of Willian & Mary. Applied & Numerical Mathematics [Managing 
Uncertainties in Multidisciplinary Optimization] 

Semyon V. Tsynkov - Ph.D., Computational Mathematics, Keldysh Institute for Applied Mathematics. 
Russian Academy of Sciences, 1991. Associate Professor, Department of Mathematics, North Carolina 
State University and Senior Lecturer, Department of Applied Mathematics, Tel Aviv University. Applied &: 
Numerical Mathematics [Active Shielding and Control of Noise] 

Michael Wagner - Ph.D.. Mathematical Programming, Cornell University, 2000. Assistant Professor. Depart- 
ment of Mathematics & Statistics, Old Dominion University. Applied & Numerical Mathematics [Applied 
Optimization] (August 2000 to July 2002) 

Gerald Walberg - Ph.D., Aerospace Engineering, North Carolina State University, 1974. President, Walberg 
Aerospace. Hampton, VA. Applied & Numerical Mathematics [Revolutionary Aerospace Research Concepts] 
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Mohammad Zubair - Ph.D., Computer Science, Indian Institute of Technology, Delhi, India, 1987. Professor, 
Department of Computer Science, Old Dominion University. Computer Science [Performance of Unstruc- 
tured Flow-solvers on Multi- processor Machines] 


X. GRADUATE STUDENTS 

David M. Bortz - Department of Mathematics, Center for Research in Scientific Computations, North Car- 
olina State University. (June 2001 to Present) 

Nathan L. Gibson - Department of Mathematics, Center for Research in Scientific Computations, North 
Carolina State University. (January 2002 to Present) 

Matthew L. Hausman - Department of Aerospace Engineering, University of Colorado. (January 2002 to 
Present) 

Gregory Hicks - Department of Applied Mathematics, Center for Research in Scientific Computations, North 
Carolina State University. (September 2000 to Present) 

Brahmadatt Koodallur - Department of Computer Science, Old Dominion University. (August 2000 to 
November 2001) 

Shiva Thatapelli - Department of Computer Science, Old Dominion University. (May 2001 to November 

2001) 
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